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ABSTRACT 

NGC 1407 is the central elliptical in a nearby evolved group of galaxies apparently destined to 
become a galaxy cluster core. We use the kinematics of globular clusters (GCs) to probe the dynamics 
and mass profile of the group's center, out to a radius of 60 kpc 10 galaxy effective radii) — the 
most extended data set to date around an early-type galaxy. This sample consists of 172 GC line- 
of-sight velocities, most of them newly obtained using Keck/DEIMOS, with a few additional objects 
identified as dwarf-globular transition objects or as intra-group GCs. We find weak rotation for the 
outer parts of the GC system (v/a ~ 0.2), with a rotational misalignment between the metal-poor 
and metal-rich GCs. The velocity dispersion profile declines rapidly to a radius of ~ 20 kpc, and 
then becomes flat or rising to ^ 60 kpc. There is evidence that the GC orbits have a tangential bias 
that is strongest for the metal-poor GCs — in possible contradiction to theoretical expectations. We 
construct cosmologically-motivated galaxy -(-dark halo dynamical models and infer a total mass within 
60 kpc of ~ 3 X IO^^Mq, which extrapolates to a virial mass of ~ 6 x lO-^^M© for a typical ACDM 
halo — in agreement with results from kinematics of the group galaxies. We present an independent 
C/iaradra-based analysis, whose relatively high mass at ~ 20 kpc disagrees strongly with the GC-based 
result unless the GCs are assumed to have a peculiar orbit distribution, and we therefore discuss 
more generally some comparisons between X-ray and optical results. The group's S-band mass-to- 
light ratio of ~ 800 (uncertain by a factor of ^ 2) in Solar units is extreme even for a rich galaxy 
cluster, much less a poor group — placing it among the most dark matter dominated systems in the 
universe, and also suggesting a massive reservoir of baryons lurking in an unseen phase, in addition 
to the nonbaryonic dark matter. We compare the kinematical and mass properties of the NGC 1407 
group to other nearby groups and clusters, and discuss some implications of this system for structure 
formation. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics — 
galaxies: clusters: general — globular clusters: general — X-rays: galaxies — 
cosmology: dark matter 



1. INTRODUCTION 

In our current understanding of the universe, 
stellar systems larger than a globular cluster (GC; 
baryon mass ~ IQ^Mq) are generally embedded 
in massive halos of dark matter (DM). Their sys- 
temic properties versus increasing DM halo mass 
show that there must be shifts between the various 
regimes of phys ical processes shaping the baryonic 
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between the scales of individual galaxies and clusters 
of galaxies, the mass of a halo's primary galaxy stops 
scaling with the halo mass, and reaches a universal 
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of the baryons become spread out as a hot intergalactic 
medium (IGM). The high-mass systems behave less 
like galaxies with their own enveloping halos and 
unique histories, and more like extended halos shaping 
the destinies of th eir captive galaxian members (e.g., 
iForbes et aD |2006a'). 

In collapsed galaxy clusters, the dominant processes 
affecting the larger galaxies involve cooling, heating, and 
stripping by the IGM, while some level of harassment 
and stripping also occurs by galaxy fly-bys within the 
cluster. However, it is thought that processes such as gas 
strangulation and major galaxy mergers occurred earlier, 
within the gr oups that went on to conglomerate into the 
cluster f e.g.. lBarneslll985t IZabludoff fc Mulchaev! 1 19981: 



iKodama fc Smai]ll200lHCarlberg et al.ll2001tlFuiital[20a 



iBrough et al.' '2006a', hereafter B-|-06a; 'Moran et alj 
|2007; ' Kawata fc Mulchaev 2008; Wilman et al. 200i; 
iKautsch et al.ll2008t IPatel et all 12009 ). In this picture, 
the characteristic properties of cluster galaxies (dom- 
inated by early-types) would primarily arise by "pre- 
processing" in groups. Thus by studying galaxy groups, 
one can learn not only abou t the ecosystems of most of 
the galaxies in the universe (jEkeliraLlllOOii), but also 
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about the pre-history of the rarer but interesting galaxy 
clusters. 

The nearby group surrounding the elliptical galaxy 

NGC 1407, sometiines c alled the "Eridanu s A" 

group jWi Umcr ct al.' '198^, 'Ferguson_&^aiidage' '1990|; 
iBrough et a l. 2006b.. hereafter B+06b), is a potentially 
informative transitional system. It has the X-ray and 
optical luminosities of a group (with 1-3 galaxies of 
at least L* luminosity, and highly dominated by early- 
types), but possibly the mass of a cluster. An early 
study of its constituent galaxy redshifts suggested a 
group mass -to- light ratio of T 2500 Tq (i?-band; 
lGouldlfT993h . and while subsequent studies have brought 
this estimate down to ~ 300-1100 (partially through 
longer distance estimates; Quintana et al. 1994; TuUy 
2005.: B-l-06b: .Trentham et al . 2006; Firth et al. .2006; 
Zhang et ani2007l hereafter Z-l-07). these are still high 
values more reminiscent of clusters than of groups. 

The NGC 1407 group may provide insight into the 
origins of the properties of galaxy cluster cores and of 
brightest cluster galaxies. The group is part of a larger 
structure or "supergroup" that should eventually col- 
lapse to form a cluster with NGC 1407 — and possibly its 
assoc iated group — at its core (B-|-06b; iTrentham et al.l 
|2006[) . Major mergers are prone to scramble the memo- 
ries of their progenitors' prior properties, but the future 
assembly of the NGC 1407 cluster appears to involve 
low mass ratios of ~1:5 to 1:10. Thus it is likely that 
the NGC 1407 group will remain relatively unscathed 
at the core of the resulting cluster, and we can make a 
fair comparison of its current dynamical properties with 
those of existing cluster cores (such as M87 in Virgo and 
NGC 1399 in Fornax) to see if these systems are part of 
the same evolutionary sequence. 

NGC 1407, like many bright elliptical galaxies, is en- 
dowed with a swarm of GCs visible as bright sources 
extending far out into its halo, beyond the region of 
readily observable galaxy light. GCs are invaluable as 
tracers of major episodes of s t ar formation in galaxies 

f^hitmore & Schweizer| 119951 : E arsen & Richtler' '200^; 
odic & Strader 2006). In NGC 1407, Ccnarro et al. 
)07l hereafter C-l-07) have used them to study the 
galaxy's early metal enrichment history, and to search for 
signs of recent gas-rich mergers. Here our focus is on us- 
ing GCs as bountiful kinematical tracers, useful for prob- 
ing both the dynamical properties of their host galaxy's 
outer parts, and the mass distribution of the surrounding 
group. 

The use of GCs as ma ss tr acers is well establish ed: see 
iBrodie fc Straded (j2006>) and lRomanowskvl ()2009l ) for re- 
views, and §5.2l of this paper for specific examples. GCs 
are so far less widely exploited as orbit tracers — using the 
internal dynamics of a globular cluster system (GCS) to 
infer its formation history and its connections with other 
galaxy subcomponents. Studies of GCS rotation have im- 
plicitly treated the GCs as pr oxies for the halo field stars 
or for the DM particles (e.g., ^Kisslcr-Patig fc Gc bhard'3 
[l998t) . but the connections between t hese different enti- 
ties have not yet been established. iBekki et al.l (|2005l 
hereafter B-l-05) presented pioneering work in this area, 
deriving predictions for GCS, stellar, and DM kinemat- 
ics from galaxy merger simulations. We will make some 
comparisons of these predictions to the data. 

The 10-meter Keck-I telescope with the LRIS multi- 



slit spectrograph (|Oke et al.lll995f ) has produced a legacy 
of the first large samples of ext ragalactic GC kinemat- 
ics outside of the Local Grou p (jCohen fc Rvzhovlll997t 
lCohen|[2000t iCote et al.l[2003l hereafter C-K 33). Now we 
prese nt the first use of DEIMOS on Keck-II f Faber et all 
l2003f l to study extragalactic GCS kinematics, represent- 
ing a significant advance in this field. Its 80 arcmin^ 
field-of-view is one of the most expansive available on 
a large telescope, and its stability, red sensitivity, and 
medium resolution make it highly efficient for acquir- 
ing accurate velocities for GCs (which are intrinsically 
red objects) even in gray skies, using the Ca II absorp- 
tion line triplet (CaT). Our results here inaugurate the 
SAGES Legacy Unifying Globulars and Galaxies Survey 
of the detailed photometric and spectroscopic properties 
of 25 representative galaxies and their GCSs. 

NGC 1407 is a moderately-rotating EO galaxy, nei- 
ther boxy nor disky, with a "core-like" central lumi- 
nosity profile, a weak central AGN, and absolute mag- 
nitudes of Mb = -21.0 and Mk = -24.8 (B-f06b). 
Its GCS displa ys a typical color b imodality of " blue" 
and "red" GCs (iPerrett et al.l ll997t [Forbes et al. 2 006U 
iHarris et all l2006l . hereafter H-^06a: iHarris ,2009l) . Its 
overall specific frequency of GCs per unit luminosity 
was previously found to be a modest Sn ~ 4, but a 
new deep wide-field study has produced a richer esti- 
mate of Sn ^ 6 (L. Spitler et al., in preparation, here- 
after S-t-09). NGC 1407 hosts bright, extended X-ray 
emission that extends as far out as 80 kpc into the sur- 
rounding group (B-f 06b; Z-l-07), although there is noth- 
ing about its X-ray pro perties that marks it as a p artic- 
ularly unusual galaxy (jOsmond fc PonmanI |200S ^ here- 
after OP04). We adopt a galaxy stellar effective radius 
of i?off= 57" = 5.8 kpc (see M.l.ll). with a distance to 
the galaxy and group of 20.9 Mpc |Forbes "eraIll2006bD — 
discussing later the effects of the distance uncertainty. 

We present the spectroscopic GC observations in ij2l 
and convert these into kinematical parameters in ^ In 
ijH we analyze the group mass profile using GC dynam- 
ics and X-ray constraints, and review independent mass 
constraints more generally. We place our findings for 
NGC 1407 into the wider context of galaxy groups in ij5l 
and Sj6] summarizes the results. 

2. OBSERVATIONS AND DATA REDUCTION 

Here we present the spectroscopic observations used 
in this paper, including new data from Keck/DEIMOS, 
and pubhshed data from Keck/LRIS. The DEIMOS ob- 
servations are outlined in and the data reduction 
in g221 The combined DEIMOS-hLRIS sample is de- 
scribed in §2.31 We also make use of i ntegrated-light 
stellar kinematics data for NGC 1407 from lSpolaor et aTl 
lj2008al hereafter S-l-08) and photometric data of the 
galaxy and its GCS from S+08 {HST/ ACS BI imag- 
ing) and_^ran_Sjj^9(^ubaru/Suprime-Cam g'r'i' imag- 
ing; iMivazaki et"al]|2002f) . 

2.1. Keck/DEIMOS observations 

The DEIMOS observations were carried out as part of 
Keck program 2006B_U024D (PI: J. P. Brodie). Candi- 
date GCs were selected from Subaru/Suprime-Cam im- 
ages of NGC 1407 taken in g'r'i' bands in very good 
seeing conditions (0.5"-0.6"). The NGC 1407 GC cata- 
log is detailed in S+09 and was supplemented with the 




Fig. 1. — Overview of spectroscopic observations around NGC 1407. Filled symbols show results from DEIMOS, and open symbols from 
LRIS. Orange circles show positions of kinematically confirmed GCs, and green squares show dwarf-globular transition objects. Purple 
triangles show borderline velocity results, red X symbols mark background galaxies, and Galactic stars are shown by blue starlike symbols. 
Small gray points show the underlying distribution of GC candidates from Suprime-Cam, down to the same limiting magnitude as the 
spectroscopically-confirmed GCs; the irregular distribution of spectroscopic GCs is caused by mask and slit placement and geometry, and 
by poor-quality regions in the imaging. The central round contour shows the fij = 23.2 mag arcsec"^ stellar isophote of NGC 1407, 
corresponding to ~ 3 R^s (2.85'). The lower-right contour shows a similar isophote for NGC 1400. The dotted line shows the approximate 
boundary of predominance between the GCSs of the two galaxies. 

target slits were aligned close to the parallactic angle, 
and had lengths of at least 5" each. Three masks pro- 
vided coverage of the NGC 1407 GC system out to large 
radii, and were positioned to avoid potential contamina- 
tion from the nearby galaxy NGC 1400 to the southwest. 
The GCS is fairly well sampled out to ^ 5', with an 
extension of ~ 10' to the northeast (see Fig.[T]). 

Spectroscopic observations were made with Keck- 
2/DEIMOS on 2006 Nov 19-21, using the "standard" 
DEIMOS setup with the 1200 1/mm grating, centered on 
~ 7500 A. This central wavelength allowed coverage of 
the CaT at ~ 8500-8700 A, independent of slit position 
for GCs at the highest expected line-of-sight velocities in 
the NGC 1407 GC system (~ 2500 km s"!), while cap- 
turing Ha for a reasonable fraction of the GCs observed. 
For this setup, our slit width of 1" gives a FWHM of 



NGC 1407 HST/ ACS GC list of lForbes et ^ ((2006bl) . 
An astrometric solution was derived for the Subaru g'- 
band image using u nsaturated USNO-B stars in the field 
(|Monet et al.l 120031 ). Objects detected on the Subaru 
imaging were tagged as GC candidates if they met color 
limits of 0.4 < g' - i' < 1.4 and 0.3 < g' - r' < 0.9, and 
if the color values were jointly consistent with the clear 
linear sequence of GCs in g' — i' vs. g' — r' color space. 

For the multi-slit mask design, GC candidates were 
ranked according to their expected S/N . No cuts were 
made on GC color beyond the initial selection of can- 
didates. DEIMOS alignment and guide stars were also 
taken from the USNO-B catalog, with a preference for 
stars in the recommended magnitude range of 15 < i? < 
20. To ensure a good alignment, at least two stars were 
placed at each end of a given DEIMOS mask. Science 
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Fig. 2.— Extracted DEIMOS spectra of NGC 1407 GCs, shifted 
to zero redshift, with the CaT region shown, and with 5-pixel box- 
car smoothing applied for visibility. The top spectrum shows a 
best-case example (#136205), the middle one is typical (#25581), 
and the lower one is a minimally acceptable example (#192009). 
The rest-frame wavelengths of the visible Ca II absorption lines are 
8498 A, 8542 A, and 8662 A. 

- 1.5A 50 km s^^ at the CaT). The three masks 
were observed for 1.5-2 hours each in individual expo- 
sures of 30 min. The seeing on aU three nights ranged 
from 0.55" to 0.8". The large number of high-quality 
spectra that will we extract in ii2.2\ is testimony that 
the ~ 16' X 5' field of view of DEIMOS, along with its 
Flexure Compensation System, provide substantial im- 
provements over Keck-l/LRIS (~ 8' x 6' field, not all of 
it usable for full spectral coverage). 

2.2. Keck/DEIMOS data reduction 

The spectra were reduced using the IDL spec2d 
pipeline written for the DEEP2 galaxy survey (J. A. 
Newman et al., in preparation). Images were flatfielded 
using an internal quartz lamp, and wavelength calibrated 
using ArKrNeXe arc lamps. After sky subtraction, indi- 
vidual one-dimensional spectra were optimally extracted. 
Some sample extracted spectra are shown in Fig. [5J 

Radial velocities were derived through cross- 
correlation around the CaT region with a library 
of standard stars spanning a wide range of spectral 
types (mid F to late K) taken both on this run and 
during previous DEIMOS runs with the same setup. 
The derived velocities varied with spectral template 
choice by only ~l-2 km s~^. In all cases the best-fit 
line-of-sight velocity was verified by visual inspection of 
both the CaT region and Ha (where possible). Small 
(< 5 km s~^) shifts in the zero-point velocities of the 
three masks were corrected using the location of the 
telluric A band. 

The velocity uncertainties were estimated from the 
width of the cross-correlation peak, and most fall within 
the range ~ 5-15 km s~^, depending on the brightness 
of the GC. To assess the accuracy of these uncertainties. 
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Fig. 3. — Comparisons of duplicate velocity measurements, 
as a function of ig-band magnitude. The black filled circles show 
repeated measurements using DEIMOS. The gray open circles show 
duplicate measurements between LRIS and DEIMOS, where we 
have corrected for a mean offset of 11 km s~^. Error bars show the 
total estimated uncertainties. In the case of LRIS, the solid bars 
show statistical uncertainties, and dotted bars include systematic 
uncertainties. 

we deliberately included duplicate observations of GCs 
on different masks. In Fig. [31 we show the results of the 
10 repeated velocity measurements (9 GCs and 1 star). 
The velocities are nicely consistent within their stated 
uncertainties, to our faintest magnitude limit of if, 23. 

Taking duplicates into account, we obtain velocities 
of 170 objects using DEIMOS, and list their properties 
in Table [TJ We include two objects serendipitously ob- 
tained from mask alignment holes, and two faint blue 
objects that are borderline cases in our conservative spec- 
tral selection criteria. These latter two 131881 and # 
160394) will be tentatively incorporated in our analyses, 
while checking that their inclusion does not affect any 
results. For a further 40 objects — generally the faintest 
in our masks — we could not obtain accurate line-of-sight 
velocities. Many of these showed peaks in the cross- 
correlation consistent with velocities expected for NGC 
1407 GCs, but because of their low S/N no clear visual 
confirmation of the CaT could be ascertained. 

Of the good spectra, 8 are clearly low-redshift galaxies 
(based on emission lines or extendedness), and another 4 
are Galactic stars (see §3.11 for star/GC demarcations). 
There are 161 GC candidates confirmed to be kinemati- 
cally associated with NGC 1407, and since 7 of the stars 
and galaxies had non-GC-like colors but were observed 
for lack of another target, our effective contamination 
rate was a mere 3%. This high yield can be attributed 
to the richness of the NGC 1407 GCS, and to the excel- 
lent multi-color Suprime-Cam imaging and pre-selection. 

One potential worry when observing the CaT is the 
flurry of bright sky lines in the wavelength vicinity. How- 
ever, we have checked carefully that either Ca2 or Ca3 
is always clear of a sky line (and usually Cal also), ex- 
cept for narrow velocity ranges around —150 km s~^, 
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TABLE 1 

Photometry and spectroscopy of globular cluster candidates around NGC 1407. 



ID 


RA 


Dec 


Rp 


PA 




ig' - i')o 


(r'-i')o 


Vp 


Aj;p 


Type 




[J2000] 


[J2000] 


["] 


n 








[km s ^] 


[km s-i] 




70163 


03:40:14.257 


-18:41:32.54 


406.1 


175.1 


22.94 


0.77 


0.22 


0.1789 




gal 


71768 


03:40:18.080 


-18:41:21.86 


403.9 


167.3 


21.93 


1.08 


0.36 


1792 


9 


GC 


75581 


03:40:08.063 


-18:40:55.52 


371.4 


188.2 


21.64 


1.00 


0.31 


1750 


13 


GO 


81722 


03:40:14.309 


-18:40:10.38 


324.4 


173.7 


22.59 


0.85 


0.29 


1439 


13 


GC 


82094 


03:40:16.221 


-18:40:07.78 


326.0 


168.9 


22.38 


1.24 


0.45 


1497 


13 


GC 


82266 


03:40:15.476 


-18:40:06.01 


322.3 


170.7 


22.16 


0.82 


0.25 


1686 


17 


GC 


84789 


03:40:18.073 


-18:39:45.97 


311.1 


163.4 


22.35 


1.16 


0.39 


1788 


12 


GC 


88899 


03:40:26.549 


-18:39:13.77 


338.5 


141.8 


21.31 


1.07 


0.36 


1524 


7 


GC 



etc. 



Note. — The first column gives the Suprime-Cam identification number, followed by the LRIS identification number from C-l-07 in 
parentheses where appropriate. Tabulated properties include the projected distance from NGC 1407 {Rp), the position angle on the sky 
(North through East), the ig-band magnitude, and two colors — as well as the heliocentric line-of-sight velocity tip ± Avp (if a galaxy, the 
redshift z is listed instead). The final column indicates whether the object is classified as a star, globular cluster (GC), intra-group globular 
cluster (IGC), dwarf-globular transition object (DGTO), or background galaxy (gal); objects with "?" have uncertain identificatio ns, and 
objects in (parentheses) have marginal signal-to-noise. The photometric values have been corrected for Galactic extinction (^Schlegel et al.l 
[HID. The center of the galaxy NGC 1407 is at 03:40:11.808 -18:34:47.92 (J2000), and the systemic velocity is 1784 km s'^. This table is 
available in its entirety in the electronic version of the article. 



+225 km s'-^, and -f 1550 km s~^. The latter value may 
be of some concern to this study, as we may selectively 
lose a few faint GCs around this velocity despite the good 
sky subtraction of the DEIMOS pipeline. 

2.3. DEIMOS+LRIS combined sample 

To increase our sample size of GC velocities around 
NGC 1407, we include the 20 velocities obtained by 
C-l-07 using Keck/LRIS. These LRIS velocities also al- 
low us to make an external check on our DEIMOS ve- 
locities, using 4 GCs in common: see Fig. |3l Given the 
stated LRIS systematic uncertainties, the duplicate mea- 
surements are fully consistent, although different instru- 
ments, telescopes, and reduction procedures were used. 
In fact, the systematic uncertainties of ^ 20 km s~^ seem 
to be overly conservative, and we use only the statisti- 
cal uncertainties, increasing them by 77% and subtract- 
ing 11 km s~^ from the velocities in order to match the 
DEIMOS data. The final LRIS uncertainties are typi- 
cally '~ 7 km s"^. In conclusion, the internal and ex- 
ternal checks give us great confidence in the reliability 
of the Keck data. Velocity accuracies of 5-10 km s^^ 
are not essential for the current study of NGC 1407, but 
they do demonstrate that DEIMOS will be suitable for 
studying the dynamics of very low-mass galaxies. 

After combining duplicate measurements, we arrive 
at a sample of 190 GC candidate velocities, including 
172 GCs and 5 other objects associated with NGC 1407 
(see >j3.1l for selection). These velocities, and associated 
source properties, are presented in Table [TJ 

3. KINEMATICS RESULTS 

Having assembled a catalog of GC candidates, we select 
the probable GCs and examine their basic photometric 
and kinematical distributions in t j3.1l We study the ro- 
tation, velocity dispersion, and kurtosis properties of the 
GCS in gSH gXl and gSIl respectively. 

3.1. GC selection and velocity distribution 



0.3 I 1 1 1 1 1 ] 1 1 1 n 1 1 1 r 




Fig. 4. — Distribution of GC candidates with velocity, 
smoothed by the measurement uncertainties and an additional 
50 km for clarity. The dot-dashed line marks the systemic 
velocity of NGC 1407. The dotted line shows a Gaussian with 
cTp = 241 km s~^. 

As all of the GC candidates have colors within the 
expected range, we first select the GCs based on veloc- 
ity. This selection is trivial, since there is a large gap 
between the lowest-redshift obvious GC (1085 km s~^) 
and the highest-redshift obvious star (239 km s~^; see 
Fig. [4]). We see no objects obviously associated with 
NGC 1400 (at a distance of ^ 12' and with a systemic ve- 
locity of ^ 560 km/s), although it is in principle possible 
for some of the low-velocity NGC 1407 GCs or the "stars" 
to be extreme-velocity GCs belonging to NGC 1400. It 
is more probable that the density of NGC 1400 GCs is 
simply very low in our spectroscopically-surveyed region 
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Fig. 5. — Color-magnitude diagram of GC candidates around NGC 1407. The symbols are as in Fig.[T] except now the small gray points 
show GC candidates only within a 10' radius of NGC 1407. The typical photometric uncertainties are shown as error bars on the left. The 
dotted line shows the adopted boundary between metal-poor and metal-rich GCs: (g' — i')o = 0.98. 

(Fig. [1] shows the photometrically-estimated cross-over 
point between the GCS surface densities of NGC 1407 
and NGC 1400), and further observations will be needed 
to determine the total radial extent of this galaxy's GCS. 

We thus derive an initial sample of 177 GCs, includ- 
ing 2 cases that are borderline because of their low 
S/N. The color-magnitude diagram of the identified 
objects is shown in Fig. \5\ superimposed on the distri- 
bution of all sources in the region from the Suprime- 
Cam imaging. Our spectroscopy samples the GCs about 
three magnitudes down the luminosity function, with 
almost equal numbers coming from the blue and the 
red subpopulations. The dividing color has been set at 
ig' ~ *')o = 0.98 by inspection of the color-magnitude 
diagram. Assuming that the GCs are uniformly very 
old and that the color differences reflect metallicity dif- 
ferences (as C-l-07 confirmed for a subsample of bright 
NGC 1407 GCs), we will in the rest of this paper 
refer to the blue and red GCs as "metal-poor" and 
"metal-rich", respectively. The GC color extrema and 
red/blue boundary of {g' - i')o = (0.70,0.98, 1.30) cor- 
respond to [Fe/Hj^v^ (-2 3,-1.1,-^0.3) for a 12 Gyr age 
(|Bruzual fc CharloOl2003D . although the absolute values 
of these metallicities should not be considered certain. 

An immediately apparent feature of the spectroscopic 
GC sample is a "blue tilt" of the metal-poor peak — 
a systematic trend of color with GC luminosity that 



has recently been discover ed in ma ny bright galaxies, 
including NGC 1407 fe.g.. iStrader e t al. 200^; H+06a ; 
iMieske et all 120061: ISpitler et ain2008a; .Harrii 120091) . 
This is now one of several detections o f such a trend us- 
ing ground-based photometry (see alsol Forte et al.|[2007l : 
iWehner et a l. 2008), which argues against clai ms that 
the blue tilt is an artifa ct of HST photometry (jKundul 
[2008tlWaters et al.l[200l . 

Our data further support a convergence of the GC 
color distribution (including the metal-rich GCs) to a 
unimodal peak at high luminosities (H-|-06a). The bright 
end of the GCLF is in many galaxies known to in- 
clude extended objects such as uj Cen, Gl in M31, 
and ultra-compact dwarf galaxies (UCDs)-more gener- 
ally "dwar f-globular transition objects" (DGTOs ; e.g.. 
Hilker et al. 1999; Drinkwater et al. 2000; Hasega n et al.l 
20051: lEvstigneeva et air[2007l: IWehner fc Harrisi 12001 7 



As the nomenclature implies, some of these objects may 
have origins and properties that are distinct from "nor- 
mal" GCs; e.g., the DGTOs could be the remnants of 
dwarf galaxies that were tidally stripped down to bare 
nuclei. It has been suggested that NGC 1407 harbors a 
sizable population of extended "GCs" (H-|-06a), and it is 
important for our current study to ascertain whether our 
kinematics sample includes DGTOs, and whether their 
dynamics distinguishes them from the rest of the GCs. 
We therefore re-examine the sizes of our candidate 
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Fig. 6.— Color distribution of GCs around NGC 1407. The solid 
curve shows the 172 spectroscopically-confirmed GCs, with shaded 
region indicating the approximate 1 cr uncertainties from Poisson 
noise. The dashed curve shows the underlying Suprime-Cam dis- 
tribution of GC candidates with the same range of radius, colors, 
and magnitude, and renormalized to match the total number of 
confirmed GCs. For clarity, a smoothing kernel has been applied 
to these curves with a width of (g' — i')o = 0.03, in addition to the 
measurement errors. 

GCs in their original discovery images. The ACS data 
(near the galaxy's center) can discern the locus of GC 
sizes near ~ 3-4 pc, and we do not support the con- 
clusion of H-|-06a; their "extended" objects seem to be 
simply examples of bright GCs with normal sizes that 
are resolved at a distance of 20-30 Mpc. We can iden- 
tify only one spectroscopic object that is a clear outlier 
(#136205, with i?cff= 8 pc), which we designate as a 
DGTO. The Suprime-Cam data are much less well suited 
for measuring sizes without a very careful PSF-based 
analysis. Again, there is one clear outlier (#160583), 
which was used for DEIMOS mask alignment but turned 
out to have a velocity corresponding to the NGC 1407 
group. It is a known galaxy candid ate named AP - 
MUKS(BJ) B033807.82-183957.5 (Madd oxet aHflQOl . 
and is well-resolved in the Suprime-Cam images, with 
i?eff ~ 25 pc. Its color of {g' — i')o = 1.04 and luminosity 
of Mi'^ = -12.5 make this object likely to be a DGTO. 

The two clear DGTOs around NGC 1407 are also the 
brightest two objects in our kinematics sample. They 
have somewhat red colors, and the second-ranked ob- 
ject was studied by C-l-07, who found nothing unusual 
about its spectroscopic properties (age 11 Gyr, [Z/H]'^ 
-0.8, [a/Fe] hO.3). Another object, #134086, cor- 
responds to the "galaxy" APMUKS(BJ) B033759.53- 
184331.2, but is unresolved in the ACS images, and we 
retain it as a GC. This gives us a first-pass sample of 175 
GCs. 

The metal-poor and metal-rich GCs are sampled fairly 
equally at all radii, which may be marginally inconsistent 
with the photometric inference that the metal-poor sys- 
tem is more radially extended; the inconsistency may be 
due to some hidden bias in the spectroscopic target se- 



FlG. 7. — Velocity distribution of GC candidates with radius. 
Blue open circles are blue (metal-poor) GCs, and red filled circles 
are red (metal-rich) GCs. The two purple filled triangles are un- 
certain (metal-poor) GC detections. The three DGTOs are shown 
by green squares. The systemic velocity is marked by a dot-dashed 
line, and the ±3cr "envelope" by dotted lines. 

lection. The distribution of confirmed GC colors is com- 
pared to the underlying distribution in Fig. \6\ The dis- 
tributions are very similar, except for an excess of mildly 
metal-rich confirmed GCs that appears as a third color 
peak. This is probably an effect of the observational bias 
toward brighter objects, which show a convergence to- 
ward color unimodality (cf. Fig. [5|) and which we will 
explore further in S+09. To a good approximation, we 
can conclude that the spectroscopically confirmed GCs 
are representative of the underlying GCS, although in 
§3.21 and §3.31 we will double-check the kinematics for 
any magnitude dependence. 

The azimuthal distribution of GC measurements is sta- 
tistically uniform over the radial range of 106"-271". 
The mean and median values of the GC velocities in 
this annulus are 1769 and 1778 km s~^, respectively. 
The mean velocity of NGC 1407 from LRIS long-slit 
stellar kinematics is 1784 km s~^ (after subtracting the 
LRIS/DEIMOS systematic offset), which we adopt as the 
systemic velocity Wsys- For comparison, the NED^ veloc- 
ity is 1779 km s"^ 

The overall distribution of GC candidate velocities was 
shown in Fig. |31 Given the uncertainties, the main 
peak of GCs is consistent with a simple Gaussian pro- 
file. We next inspect the velocity distribution with radius 
(Fig. [7]), seeing no obvious asymmetries. 

There is one GC (#111003) with an "outher" veloc- 
ity (2692 km s""'^), well separated from the rest of the 
distribution. Calculating a smoothed profile of velocity 
dispersion with radius (see i )3.3p . we define a 3 cr en- 
velope outside of which there is unlikely to be any ob- 

^ http://nedwww.ipac.caltech.edu/ . The NASA/IPAC Extra- 
galactic Database (NED) is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 
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TABLE 2 

KiNEMATICAL PROPERTIES OF NGC 1407 GCS IN VARIOUS SUBSAMPLES. 



Subsample 


A'gc 


''rot 

[km s"-*-] 


do 

[°] 


Significance 


o-p 
[km s"-*-] 


fRMS 
[km s"-*-] 


f rot / CTp 




All radii 
All, 0.35'-ll.l' 
Blue, 0.35'-ll.l' 
Red, 0.35'-ll.l' 


172 
86 
86 


30+?!' 


47+5S 

— 47 

129t« 


47% 
44% 
68% 


^^'-'—12 

232l?g 
243t?i 


241+;* 

— 12 

234i?i 
2471?^ 


0.13 + °?3 
—0.13 

24+°' ^3 
'^■^^-0.18 


-0.32 ± 0.37 
-0.32 ± 0.51 
-0.28 ± 0.51 


Small radii 
All, 0.35'-1.73' 


33 


00+62 

— 80 


-178+^1 

— 48 


57% 


291+tn 

— 30 


298+t? 

— 31 


0.34+°'?2 

—0.28 


-0.61 ± 0.80 


Azimuthally complete 
All, 1.77'-4.52' 
Blue, 1.77'-4.52' 
Red, 1.77' -4.52' 


75 
38 
37 


70+29 
'^-37 
19+33 


_-,2+24 
-^^-25 
-115+^°'' 


93% 
3% 
98% 


20ll« 
198t?^ 

190l?9 


209t?° 
1981?^ 
220I22 


'^•''^-0.19 
06+°" 
n 7q+0'22 
' ^ — 0.28 


-0.38 ± 0.55 
-0.26 ± 0.75 
-0.41 ± 0.76 


Large radii 
All, 4.55'-ll.l' 
Blue, 4.55'-ll.l' 
Red, 4.55'-ll.l' 


64 
36 
28 


112+*''' 
^^^-112 
gi+40 


io2j:« 

-,,+46 
^^-54 


74% 
78% 
47% 


237+24 
257lfg 
194132 


272138 
198133 


30+°'''' 
U.^>U_(j 23 

44+° i9 
'^■^^-0.44 

32+° 2i 
"■■^^-0.32 


-0.36 ± 0.59 
-0.59 ± 0.77 
-0.24 ± 0.86 


Outside center 
All, 1.77'-11.1' 
Blue, 1.77'-11.1' 
Red, 1.77'-11.1' 


139 
74 
65 


c:9+26 
^^-30 
54+29 

11^+32 


32 ±29 

122+**^ 

9+17 
^-19 


83% 
62% 
99% 


222lj^ 
233+22 
1941?° 


22511^ 
237+?g 
21ll?i 


23+°'^2 
^•^■^-0.13 
23+°-i3 
"■^'-'-0.23 
59+° !^ 

U.OM_Q 21 


-0.34 ±0.41 
-0.32 ± 0.55 
-0.39 ± 0.59 



ject (if the distribution is Gaussian). The rogue GC 
is far outside of this envelope. It is a fairly bright, 
metal-poor object, which is borderline extended in the 
Suprime-Cam images, and whose velocity should be re- 
hable. We conjecture that it belongs to either a loose 
intra-group population of GCs rather than to the GCS 
immediately surrounding NGC 1407 (cf. the " vagrant" 
GC in the Fornax cluster; ISchuberth et al.|[2008D — or else 
to a sparse, dynamically distinct subpopulation of GCs 
that were propelled outward by three-body i nteractions, 
or were stripped from infa lling satellites fe.g.. lSales et al.l 
I2007t iLudlow e t al. 2009). This object is omitted from 
the following kinematical analyses®. 

There is a second object (=1^135251, 2505 km s~^) that 
lies just below the upper 3 a boundary. It is one of the 
brightest and reddest GC candidates around NGC 1407, 
appearing almost as an outlier in the color-magnitude di- 
agram (Fig. [5]); from the Suprime-Cam images, it may be 
unusually extended. Noticing that one of the previously- 
identified DGTOs also had an extreme (blue-shifted) ve- 
locity, we suspect that there is a sizable population of 
unrecognized DGTOs concentrated toward the galactic 
center, with a high velocity dispersion. Indeed, a dispro- 
portionate share of the borderline extended objects seem 
to be red objects with high relative velocities, residing 
at small radii. Therefore we assign this third object to 
the DGTO category, and keep in mind that the GCS 
kinematics inside ~2'~3' (~ 15 kpc) may still be "con- 
taminated" by a subpopulation of DGTOs (cf. the case 
of NGC 1399, w here a dozen DGTO s have been found 
within - 15 kpc: iThomas et al.ll2008[) . 

8 Neither this object nor any other GC candidate has a velocity 
that would place its pericenter near the virial radius as found in 
m which implies that they are all associated with the NGC 1407 
group. 



We also note that #131881 has the lowest velocity in 
the sample (1085 km s~^), but we originally classified 
it as a borderline spectrum, so we remove it from our 
sample. Our final sample of GCs thus has 172 objects. 

3.2. Rotation 

As a rough measure of the amount of rotation in the 
NGC 1407 GCS, we perform a least-squares fit to the 
velocities as a function of azimuth, fitting the expression: 

v{e) = Vsys + Wrot sin(6'o - 6*), (1) 

where v-^ot is the rotation amplitude, and 6q is the di- 
rection of the angular momentum vector. T he physical 
interpretation of this model is discussed in ICote et all 
()2001l hereafter C+01; note the different sign conven- 
tion). C+01 and some other studies have used error- 
weighted least-squares rotation fitting (equivalent to 
fitting). However, this approach is not so appropriate 
for a system with an intrinsic dispersion, and it is prefer- 
able to use straightforward least-squares fitting (equiv- 
alent to minimizing the rotation-subtracted dispersion). 
The differences become important when the true under- 
lying rotation is weak ( yrot/cp $ 0.4). Further refine- 
ments are poss i ble (e .g., iKissler-Patig k, Gebhardtl[T998l : 
iWoodlev et a^[2007^ ■ but for now we will keep things 
simple using least-squares fitting. The procedures we de- 
scribe below for evaluating statistical uncertainties and 
significance follow Sharpies et al. (199S). 

Using the full data set of 172 GCs, we find Urot = 
30I3Q km s^^ and 6*0 = 47I47 degrees, where the uncer- 
tainties have been estimated through Monte-Carlo simu- 
lated recovery of the best-fit models. An additional test 
with random reshuffling of the GC data point position 
angles shows a 53% probability of a rotation amplitude 
this high happening just by chance, i.e., there is no sig- 
nificant rotation detected. The position angle is close 
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Fig. 8. — Rotation of NGC 1407 GCS, illustrated by velocities 
plotted against position angles. The symbols are as in Fig.[3 with 
GCs in only the azimuthally- complete region (1.77'— 4.52') plotted. 
The best-fit sine curves are shown for all GCs (solid), metal-poor 
GCs (dotted), and metal-rich GCs (dashed). For illustration, the 
± 2 (Tp boundaries of the metal-rich GCs are also shown, corre- 
sponding to the dispersion of this subsample. The horizontal line 
shows the systemic velocity, and the vertical lines mark the galaxy 
photometric major axis (~ 55° at 3'; S+09). 



to the galaxy's major axis (PA 2± 55-65 degrees; S-l-08; 
S-l-09), implying that the GCS may rotate around Taiher 
than along the major axis. 

Before drawing any conclusions, we must recognize 
that the rotation field can be complicated by a mixture 
of metal-poor and metal-rich GCs with very different ro- 
tation patterns, by variations with radius, and by po- 
tential DGTO contamination in the innermost regions. 
Therefore we fit the rotation for different subsamples of 
metallicity and radius, as summarized in Table [H The 
most secure result is for the azimuthally-complete region 
of Rp = 1.77'-4.52' (11-27 kpc), where there is highly 
significant rotation close to the major axis, produced by 
the metal- rich GCs (see Fig. |8]). At larger radii, there 
is not a clear rotation signature for the metal-rich GCs, 
but the metal-poor GCs appear to rotate around an or- 
thogonal direction intermediate to the major and minor 
axes. The combination of the misaligned metal-rich and 
metal-poor rotation vectors outside the center produces 
a possible overall rotation around the major axis (with 
Vrot/cTp ~ 0.1; see §3.31 for the calculation of the velocity 
dispersion tXp). Looking in more detail at the kinemat- 
ics, the rotation of both subpopulations is driven by the 
region Rp ~ 2'~6' (13-38 kpc). 

In addition, we find evidence for a cold moving group 
of metal-poor GCs at ~ 10' 60 kpc) to the Northeast, 
with Vj-ot = 224^^2 km s~^ and CTp = 14tg° km s^^ (this 
feature is visible in Fig.[7|). These GCs have a wide range 
of luminosities but a relatively narrow spread in color 
[{g' — i')o — 0.89 ± 0.06] , suggesting a shared formational 
history. One speculative interpretation is that the cold 
group traces the recent disruption of a disk galaxy (cf. 



iFeldmann et al.ll2008f) . 

We have further investigated dependencies of rotation 
on GC luminosity, and while it at first glance appears 
that the brighter GCs rotate more strongly and with a 
different axis than the faint GCs, on closer inspection, 
this may be simply due to spatial biasing with magni- 
tude of our selected GCs. The re-introduction of the 
DGTO velocities (see would change some of the ro- 
tation vectors noticeably. Given the focus of this paper 
on the mass analysis, we will defer closer examinations 
of rotational properties to a future paper with a more 
complete data set. 

We would like to make one final comparison, between 
the rotation of the GCs and the field stars in NGC 1407. 
The general hypothesis is that the metal-poor GCs in 
galaxies are closely related to their metal-poor halo stars 
(which are observationally inaccessible at the distance of 
NGC 1407), while the metal-rich GCs are related to their 
old metal-rich bulge stars (which correspond to the bulk 
of the diffuse st arlight in an old e lliptical like NGC 1407; 
seelPritzl et al.l i2005: Forte et al.i ,2007: Baum gardt et all 
l2008f ). Unfortunately, the existing stellar kinematics 
data in NGC 1407 extend to only ~ 40" along the major 
axis (S-l-08), which is the region with only minimal num- 
bers of measured GC velocities, and furthermore with 
possible DGTO contamination. Comparing the stellar 
and GC rotation fields will thus require further efforts, 
either by extracting galaxy rotation direc tly from the 
same DEIMOS slits used for the GCs (cf. iH orriseialJ 
120081 : iProctor et aT]l2008l : iProctor et aT]l2009f ). or by ob- 
serving planeta ry nebulae (PNe) a s proxies for metal- 
rich halo stars omanows k^ l2006f) . We do note that 
the field stars and metal-rich GCs in the overlap region 
were found bv lSpolaor et al.l (|2008b( l to have similar ages, 
metallicities, and a-element enhancements. If the metal- 
rich GCs are directly associated with the field stars, then 
the increase of Wrot from ~ 50 km s~^ in the central parts 
(S-l-08) to ~ 100 km s~^ in the outer parts may support 
the hypothesis that th e "missing " angular momenta of 
elliptical galaxies (e.g., 'FaU"l983) can be found in their 
halos (e.g.. Weil fc Hcrnquist 1996; B+05). 

3.3. Velocity dispersion 

We construct a projected GCS velocity dispersion pro- 
file from these data out to a distance of 60 kpc from 
the central galaxy. Technically, we compute the RMS 
velocity profile, wrms ^ ^jYji{vi — Vgysy/N, which in- 
cludes the (small) contribution from rotation, but we 
will refer to it freely as the "velocity dispersion" . The 
true dispersion cTp is also calculated by subtracting the 
contribution from the rotation f tj3.2p . and both wrms 
and CTp are shown in Table [21 We compute the binned 
disp ersions and its unc ertainties using standard formu- 
las (|D anese et al.lll980( ). and complementarily construct 
smoothed dispersion radial profiles using maximum like- 
lihood fitting to a Gaussian line-of-sight velocity distri- 
bution, weighted at each point by neighboring data (see 
iBergond et al.|[2006f ). We also use this maximum likeli- 
hood technique to fit power-law functions to the disper- 
sion profiles: 

4ms(^p) = 4ms,o X (Rp/Roy-. (2) 

With the full data-set of 172 GCs, there is a strik- 
ing transition in the dispersion profile (Fig. [9]) from a 
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Fig. 9.— Projected RMS velocity profile of NGC 1407. The cir- 
cles with error bars show binned GCS values with no color selection. 
The solid curve shows a smoo thed GCS p rofile using a variable- 
width Gaussian kernel (see iBergond et al ll2006i for details). The 
shaded region indicates the l-cr uncertainty boundaries for the full 
GCS sample, as computed by Monte Carlo simulations. The dotted 
and dashed lines show the profiles for the metal-poor and metal- 
rich GC subsamples, respectively. The star symbols with error 
b ars show long-s lit stellar kinematics data from S-l-08, as updated 
in lProctor etaLl (2009). 

rapid decline inside 1.8' (11 kpc — confirming the finding 
of C-f 07) to a fairly constant value outside of this radius 
(7p — —0.07 ± 0.26). There is also the hint of another 
decline outside of 7.5' (50 kpc), but it is not significant. 
The constancy of the dispersion profile out to ~ 10 Rcs 
is immediately a strong indication of a massive DM halo 
around NGC 1407. 

To better understand the dispersion profile, we analyze 
the metal- rich and metal-poor GCs separately (Fig. ^ . 
The dispersions of the metal-rich and metal-poor clus- 
ters are similar over the range of Rp ~ 3'-4.5' (~ 20-30 
kpc). Inside this range, the metal-rich GCs have a signifi- 
cantly higher dispersion than the metal-poor GCs, and at 
larger radii, their dispersion is lower. This difference was 
already visible by eye in Fig. [7l and translates into very 
different slopes with radius, with 7p = —0.26 ± 0.23 {in- 
creasing) for the metal-poor GCs, and 7p — -1-0.54 ±0.19 
(decreasing) for the metal-rich GCs. 

Such a difference in dispersion slope is qualitatively ex- 
pected from the differences in spatial density profiles of 
the two subpopulations (the blue GCS is more extended 
than the red GCS), but explaining this situation quanti- 
tatively may be a challenge (investigated in i j4.1.4p . The 
steepness of the inner slope (driven by the metal-rich 
GCs) would have been even steeper if the probable DG- 
TOs had been left in the sample (see WS.l^i . which we 
surmise may still include contamination from a popula- 
tion of as yet unresolved, metal-rich DGTOs on fast or- 
bits at small galactocentric radii (we will return to this 
issue in H5.l^ . Within the uncertainties of our limited 
data set, we do not find any trend of dispersion with GC 
luminosity. 



Since we are unsure about the nature of this steep in- 
ner slope, we will omit the inner dispersion data point 
from our mass profile analyses ( ij4.1l) . assuming that it 
is affected by a distinct subpopulation of objects with 
very different orbits than the rest of the GCs. Note that 
the dispersion of the metal-rich GCs is inconsistent with 
the field stars' dispersion in the narrow region of overlap 
(Fig.®. 

3.4. Velocity kurtosis 

Other measures of line-of-sight velocity distribution 
(LOSVD) beyond rotation and dispersion can provide 
key information about the orbital structure. In par- 
ticular, measures of departures of the LOSVD from a 
Gaussian shape, such as provided by the fourth-order 
velocity moment, are important diagnostics that can be 
crudely used to distinguish between radial and tangen- 
tial orbits. Here we compute a basic projected reduced 
kurtosis statistic that is approximated by: 



1 / 
N ^ 



^RMS 



3± 



(3) 



with a more complicated bias-corrected expression for 
small- number statistics given by G2 in iJoanes fc Gil] 
(|1998l )Q. For a Gaussian LOSVD, Kp = 0, which de- 
scribes isotropic orbits in a logarithmic potential. The 
kurtosis is only useful for regions over which v-^ms is ap- 
proximately constant, so we exclude the inner regions of 
the GCS {Rp < 1.73') from our analysis. 

For the overall GCS outside the central regions, we find 
Kp = —0.30 ±0.41, i.e. a suggestion of negative kurtosis. 
The "fiat-topped" behavior of the platykurtic LOSVD 
can be marginally seen by eye in Figs. [4] and [3 where 
the velocities do not peak strongly around Wgys (partic- 
ularly at large radii). The LOSVD is still statistically 
consistent with being Gaussian (the kurtosis estimator 
is unfortunately not a very powerful measure of LOSVD 
shape), but a strongly-peaked (leptokurtic) distribution 
does appear to be ruled out, with important implications 
for the orbits as we will see in the next section. For the 
metal-poor GCs, we find Kp = —0.24 ± 0.55, and for the 
metal-rich GCs, we find Kp = —0.30 ±0.59, so there is no 
overall difference detected in the LOSVD shapes between 
the two subsystems. 

We would like to further consider the behavior with 
radius Kp(i?p) but the data are inadequate for detecting 
significant trends. We do see signs of the kurtosis be- 
coming stronger with GC luminosity: see Fig. IIOI where 
there appears to be a dearth of velocities near Vgys for 
i'g <, 21.3. This effect is driven by bright, central ob- 
jects, and does not impact the result of negative kurtosis 
outside the central regions. 

One caveat about our entire kurtosis analysis is that a 
Gaussian LOSVD was assumed in i )3.1l to remove outliers; 
for example, restoring #111003 to the data-set would 
imply a higher vbms and a significantly positive Kp. More 
robust measures of LOSVD shapes and outlier removal 

" In principle, Kp can be compromised by velocity measurement 
uncertainties, rotation, and small-number statistics. However, we 
have modeled these effects analytically and with Monte Carlo sim- 
ulations, and determined that they are negligible for this data-set, 
except for the small number statistics which may bias the results 
by ~ -0.1. 



Mapping the Dark Side with DEIMOS 



11 



2500 



2000 



1500 



1000 



• o o » 

._._.o._..°,...._^ ° 



19 



20 



21 



22 



23 



Fig. 10. — Distribution of GC candidate velocities versus mag- 
nitude. Symbols are as in Fig. [7] Note the paucity of bright GCs 
near the systemic velocity. 

will eventually require detailed self-consistent dynamical 
modeling. We will return to the dynamical implications 
of the kurtosis in ij4.1.41 and to the issue of the luminosity 
dependence in 

We also note that there are no faint GCs {i'g > 22.5) 
with V ~ 1550 km s~^, which may be due to systematic 
loss near this velocity from sky line collisions fi i2.2p . But 
even if this is the case, we would expect to have lost only 
about one GC to the effect, which should not compromise 
our results. Additionally, a few GC candidates near the 
extreme mask edges experienced wavelength-dependent 
vignetting that was not recognized during the original 
mask design. However, we have verified that this issue 
would not have caused any velocity bias in the recovered 
GC spectra. 

4. NGC 1407 MASS PROFILE 

Reaching the main focus of our paper, we consider var- 
ious aspects of the mass profile of the NGC 1407 system. 
In iJ4.11 we model the dynamics of the GCS to derive a 
set of permitted mass models. In tj4.21 we present an 
X-ray-based mass analysis, and in ^4.31 we compare in- 
dependent mass constraints in NGC 1407 and in other 
galaxies. 

4.1. GCS dynamics 

Here we model the dynamics of the GCs around 
NGC 1407, beginning with a summary of the model- 
ing assumptions in i j4.1.11 then results with no DM in 
HAA/Zl and a ACDM-motivated model in H4.1.'6\ Alterna- 
tive models that better reproduce the data are explored 

in gm 

4.1.1. Modeling assumptions 

Our approach is to construct a sequence of 
cosmologically-motivated galaxy-|-halo mass models, cal- 
culate their projected GCS RMS velocity profiles 



WRMs(-Rp)i ^iid compare these predictions to the data. 
For this initial dynamical model of the GCS of NGC 
1407, we make several simplifying assumptions. First 
is that the system is spherically symmetric. The 
NGC 1407 outer galaxy isophotes have an ellipticity of 
e ~ 0.05 (S-l-08), supporting sphericity as a reasonable 
approximation — which seems to be the case in general for 
massive ellipticals at group and cluster centers, although 
flattening may become importan t at radii of ^ 20 kpc 
(|Portcr et al.l il991t [Rvdenl ll996l. 

Secondly we assume a particular form for the GCS 
velocity dispersion tensor that is nearly isotropic near 
the central galaxy, and becomes increasingly radially- 
anisotropic at larger galactocentric radii (details below) . 
We also assume that the system is in dynamical equilib- 
rium and that the solutions to the Jeans equations are 
physical and stable (these assumptions are supported by 
analyse s of simulated DM h alos with no recent major 
merger; iDiemand et"am2005l hereafter D+05). 

Thirdly we represent the cumulative mass profile MM 
by a constant mass-to-light ratio (T*) Sersic, (19681) 
model for the galaxy's field stars, plus a DM halo char- 
acterized by an "NFW" profile, expressed as a circular 
velocity that is independent of distance and bandpass: 



2, ,_GM(r) 



1 - 



r[(3-p)m, (r/as)V"'] 



T[{3-p)m] 



r 



1 



1^ / ' (4) 
1 + r/rs 

where m and flg are th e Sersic index and scale r adius, and 
p is a function of m (|Mamon fc Lokasl[2005bl . hereafter 
ML05), while Vs and are the characteristic velocity 
and spatial scales for the DM halo (N avarro et al.i,199&) . 
The Vc profile of the halo reaches a maximum of 0.4:65vs 
at a radius of 2.16rs. 

Our Sersic stellar mass model is obtained by fitting a 
preliminary combined photometric profile (/-band ACS 
-|- i'-band Suprime-Cam) outside of a radius of 1", re- 
sulting in fit parameters corresponding to (m = 4.32, 
i?eff= 57" = 5.8 kpc. Mi = -23.3). As is typical for 
giant ellipticals, these values are fairly uncertain (total 
luminosity by at least 20%, and i?cff by at least 30%) 
because of the large fraction of light contained in the 
faint stellar halo. One obtains a f ar superior fit to the 
galaxy light using a Nuker model (|Lauer eTalllTMl or 
a core-l-Sersic model (e.g., iTrujillo et al.ll200'^ . but the 
former does not converge to a finite luminosity, and the 
latter is not simple to incorporate in the modeling equa- 
tions. At the radii corresponding to our GC velocity mea- 
surements (to ~ 10'), the total luminosity of the Sersic 
model deviates from that of the core-|-Sersic model by up 
to 15%, but as we will see, the total mass in this region 
is dominated by DM, and in any case, the luminosity 
profile is not well constrained outside of 2'. 

The galaxy mass profile is uniquely determined by 
the stellar mass-to-light ratio T*, which we estimate 
from a stellar populations analysis (Z-l-07) to be T* j = 
2.12 (T*,s — 4.45 Tq^b), assuming a Kroupa 

IMF^o. The adopted NGC 1407 galaxy mass is thus 

^'^ This value for T* may seem low relative to historical findings 
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Fig. 11.— NGC 1407 GCS surface density profile. Points with 
error bars show Suprime-Cam data with a background contam- 
ination level subtracted (S+09). The solid curve is our best-fit 
Sersic model. For comparison, the NGC 1407 galaxy light model 
is shown as a blue dot-dashed line (with arbitrary normalization). 
The other curves show profiles corresponding to the ACDM simula- 
tions of D-l-05, for overdensities of 2.5 a (dashed) and 4 cr (dotted). 



= 2.0 X 10^^ Mq] we will examine the impact of T* 
on our results in §4.1.41 and i i4.3.1l 

The density profile for the GCS is obtained by a Sersic 
model fit to preliminary surface density data Sgc(^p) 
from Suprime-Cam, over the radial range of 0.8'-23' 
(S-f09). As shown in Fig. [11] the fit is characterized by 
the parameters (m = 1.39 ± 0.09, R^s^ 3.52' ± 0.07' = 
21 kpc), demonstrating that the radial extent of the GCS 
of NGC 1407 is far greater than that of its field starlight. 
Dividing the GCS into subpopulations, the metal-poor 
component is more extended than the metal-rich, with 
(m ~ 1.3, i?cff ~ 4.3') and (m - 1.2, i?efr ~ 2.8'), re- 
spectively. We have also verified that these density fits 
are fair representations of the bright GC subpopulation 
probed in our kinematics sample (jq < 23). The uncer- 
tainty in the GCS density profile is not explicitly included 
in our dynamical modeling, but we will explore the effects 
of reasonable variations in this profile as needed. 

The model GCS dispersion profile urms(^p) is calcu- 
lated with the spherical Jeans equations as summarized 
in ML05^^. In using these equations, we are also as- 



in oth er bright ellipticals [e.g., T, / ~ 3 Tq,/ in |Cappella ri et al.l 
I^OOB ). and T«,b ~ 6 Tq^b in [Gerhard et al. (2001)], and 
for NGC 1407 in part icular, the stellar populations models of 
IHumphrev et al.1 (12003 . hereafter H-|-06b) imply T^^b = 5.7 ± 
0.5 Tq B. However, the Z+07 value is meant to characterize the 
stars outside of 0.5 R^ff rather than in the oft-studied central re- 
gions (r^ 0.1 -Rcff in the case of H+06b); ano ther study ave raging 
over 1 Roff found T*_s ~ 4 for NGC 1407 (ITortora et al.ll2009il . 
The full range of possible values Z+07 found for different radii 
and modeling assumptions was T* s = 3.8-7.1 s for a Kroupa 
IMF, and up to 1 1.2 Tq,b for a S'alpeter IM F. 

We use the 'Prugniel fc SimienI ([1997') function p(m) in the 
GCS density profile since i n this case it appears to provide a bet- 
ter match than the iLima Neto et al] ([l993) function to the Sersic 
surface density profile. 




/?p [arcmin] 

Fig. 12.— NGC 1407 GCS projected RMS velocity profile. 
Points with error bars show the measurements, while curves show 
model predictions that are described in more detail in the main 
text. The black dotted curve (S) shows a stars-only galaxy (no 
DM). The green short-dashed curve (GR) is a model with a ACDM 
halo and a theoretically-motivated GC anisotropy profile (Eq. (B]!. 
The dark blue dot-short-dashed curve (GI) is a similar model with 
an isotropic GC system. The orange solid curve (II) is a power-law 
fit to the data. The red solid curve (GT) is a ACDM halo with 
tangential anisotropy. The purple long-dashed curves come from 
the X-ray based mass model: with either /3 = —0.5 (XT) or an 
anisotropy that varies from /3 = 1 in the center to /3 = — oo at 
large radii (XG; see Fig. I13I I. 

suming that the net rotation can be simply folded into 
'^RMS along with the velocity dispersion, since the ob- 
served rotation in NGC 1407 is not dynamically impor- 
tant (wrot/fp ~ 0.4: see H3.2^ . We generally attempt to 
fit the dispersion data outside of 1.8' only, since the cen- 
tral values may be skewed by DGTO contamination (see 
i j3.3[) . We will return to this assumption in i j4.1.4l 

4.1.2. Constant mass-to-light ratio models 

Our first model is a galaxy with no DM halo, or a 
"constant-T" model; see Table [3| for the details of the 
main models that we will discuss. In this case we keep 
things simple by adopting an isotropic velocity dispersion 
tensor: /3(r) = 0, where f3 = 1 — ag/a^ specifies the bias 
of specific kinetic energy in the tangential direction ag 
relative to the radial direction a^. This model's predicted 
WRMs(-Rp) for the GCs in the outer parts is much lower 
than the observed values [see model S in Fig. [T^l with 
comparison data corresponding to the v^MsiRp) profile 
from i i3.3j . In principle, the data could be roughly fitted 
by setting the distance to 40 Mpc, the stellar mass-to- 
hght ratio to T*j ^ 7 Tqj (T*,b ~ 15 Tq^b), and the 
outer GC orbits to nearly circular (/3 ~ —10), but this 
solution is highly contrived and implausible. 

The discrepancy between the velocities measured and 
predicted for a simple stellar galaxy is strong evi- 
dence for either non-classical gravity beha vior (e.g., 
iBrownstein fc Moffatl[200l [Angus et al.ll2008D . or for an 
unseen mass component — or for both. It is beyond the 
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TABLE 3 

Summary of noteworthy mass models for the NGC 1407 group. 



Label 




/3{r) 




Pa 


rs 


^ vir 




A/vir 


Xgc 




Model description 




L W^-cs J 




[^] 




[kpc] 


[Mpc] 




flO^^Afol 








S 


4.45 











0.18 




0.03 


43.5 (4) 


46.0 (7) 


Constant-T 


GR 


4.45 


Eq. [6] 


imotTo 0.83tlil 230l??o 


-, oc+0.65 
-"-■•^•^-O.ST 


c Q+2.7 
^•^-2.0 


1 , 0+31.9 


2.6 (2) 


10.9 (7) 


Radial AGUM 


GI 


4.45 





-,000+600 


1 2+l'^ 


170+2"" 


-,0+0.52 
^•^^■^-0.28 


„ p,+3 1 
6-7I2.3 


„ .+17 7 

8-4i5:o 


1.1 (2) 


11.1 (7) 


Isotropic ACDM 


II 



















0.8 (2) 


81.7 (7) 


Isotropic GC inversion 


GT 


4.45 


-0.5 


iZUU_2oo 


1 5+2'^ 


140l^g» 


1 oo+°--*5 

i.UU_Q 24 


7 0+3.4 
'■•^-2.5 




0.8 (2) 


11.6 (7) 


Tangential ACDM 


XR 


4.45 


Eq.[6] 


1300 ± 100 


11+4 


52+1^ 


U-«»_o.06 


17 ±2 




8.4 (4) 


0.8 (5) 


X-ray based NEW 


XT 




-0.5 








^•'°-0.04 






5.6 (4) 


0.0 (2) 


Tangential X-ray based 


XG 




Fig. [13] 








'"-0.04 






0.7 (4) 


0.0 (2) 


GC/X-ray consonance 


C 


4.45 


-0.5 


1300 ± 200 


2 4+1.8 


120 ± 40 


1 07+°-^** 
-^•'^'-0.17 


9 1+2-3 




4.0 (2) 


7.2 (5) 


Gonsensus ACDM 


GG 


4.45 


-0.5 


1400+2°" 71+060 OSO+I"" 


1.26 ± 0.21 






0.7 (2) 


14.0 (7) 


GCs + galaxies 



= v'^ / {4:nGr'^). The reduced statistic, = X^/^^doit is provided along with 



Note. — The characteristic NEW halo density is ps 
the number of degrees of freedom Ndof- 

scope of this paper to treat alternative gravities, so we 
will simply presume to be witnessing the effects of a mas- 
sive DM halo. Note that the missing mass is not at- 
tributable to the hot gas in the system, whose mass of 
^ lOi"Af0 (Z+07) is less important than even the stellar 
component. 

4.1.3. Fiducial ACDM models with anisotropy 

We next examine a fiducial set of models including DM 
halos. These models are parameterized by the virial mass 
Mvir and the halo concentration Cvir = r^h/Ts, which 
are correlated according to the predictions of A^-body 
cosmological simulations by the relation: 

Cvir = 17.88 X (Mvir/10"Mo)-"-i25, (5) 

with a scatter in Cvir (at fixed Mvir) of 0.14 dex (though 
the scatter ma y be smaller for high-mass, relaxed ha- 
los: iBullock et al 2001: W echsler et al.ll2002l : see §3.2 in 
iNapolitano et al]|2005l ). Here we have adopted a virial 
overdensity of Avir = 101, a power-spectrum normaliza- 
tion of (jg = 0.9, and a Hubble constant of — 70 
km s~i Mpc~i. T he relation ^ was derived from 
IBullock et all (|2001f ) for galaxy -mass halos, but turns out 
to also be valid for higher mass es, differing by at most 
5% from the updated relation of iNeto et al.l (|2007| ) over 
the range Mvir = (2-70) X lOi^Mg. Lower values for 
(Tg would imply systematically lower concentrations, but 
the effect is probably smaller than th e intrinsic scatter 
of the concentrations (jDuffv et al.ll2008( ). 

For simplicity, we do not include any back-reaction of 
the baryon s on the DM profile (such as adiabatic contrac - 
tion: e.g.. iBlumenthal et al.iri986t iGnedin et al.ll2004f ). 
which would be only a small correction at the outer DM- 
dominated radii that we are probing (see also H-f06b). 
We also do not include the most recent refinements to 
the NFW profile used in equation |4] (e.g., iNavarro et al.l 
[200llGaoet"aI] [20081) . 

The anisotropy profile is motivated by numerous the- 
oretical studies that generically predict that DM ha- 
los, and the baryonic tracers within the m, are radi- 
ally anisotropic in their outer regions (e.g., Ivan Albadal 



1981 [Cole^ fc Lace"v^ IT9961 iDiaferid 119 91 ^^'et^ 
20001: iFuku shigc & Makino 2001; Vespcr ini etlll 1 20031: 
Rasia et al.ii2q04i: iWoita k et al. 2005; Dckel et al.''200a 
ML05; D+05; lAbadi_et al. 2006; Hansen & Moore 200a 
[Sommer-Larsen 2006; Nipoti et al. 2006; Onorbe et al] 
1200 7; McMill an ct al. ,2007; Prieto fc Gnedin 2008) . 
D-l-05 in particular explored the dependence of the 
anisotropy of tracers on their concentrations within their 
host halos, and we adopt their equation (8) to represent 
the anisotropy profile of the NGC 1407 GCS^^: 

/3(r)=/3o^— . (6) 

To assign values to the asymptotic anisotropy /3o and 
the break radius tq, we first compare our GCS surface 
density profile to the parameterized projected profiles of 
the D-l-05 simulations corresponding to DM subsets char- 
acterized by an overdensity of n times the RMS back- 
ground fluctuation (see their equation 1). We ignore the 
central cusp in these density profiles since D-l-05 did not 
report the natural variation around a typical log slope 
of —1.2, and since the central GCS number density may 
be altered considerably by destruction processes. Focus- 
ing on the outer density profiles, we find that n ~ 4 is 
needed to reproduce the steep decline seen in NGC 1407 
(see Fig. ITU) . The data include both metal-rich and 
metal-poor GCs, so this fit is not appropriate for de- 
ter mining GC fo r matio n redshift as proposed by D-l-05 
and iMoore et al.l (|2006f ). but rather provides a heuristic 
reasonable assumption for the anisotropy profile. The 
log density slope at rg of the n = 4 simulation model is 
—3.3, and finding the location in our GCS density model 
with the same slope, we arrive at ro — 6.33' = 38.5 kpc. 
As specified by D+05 equation 8, n = 4 corresponds to 
/3o = 0.8; Fig. [T3| shows the resulting anisotropy profile. 
The values of /3o and tq are not very sensitive to our 
choice for n. 

^'^ The D+05 anisotropy profile has the rotational component 
subtracted, unlike our modeling convention, but this difference im- 
pacts their predictions at the level of A/3 ~ —0.02 at most. 
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Fig. 13. — Velocity dispersion anisotropy profiles for GCs 
in NGC 1407. The black horizontal dotted line shows isotropy. 
The green short-dashed line shows a theoretically-motivated pro- 
file (Eq. [GJi. The orange horizontal solid line shows the dynamics- 
based estimate of the GC anisotropy ( ^14. 1.41 1. The purple long- 
dashed line with surrounded shaded region sho ws the profile nec- 
essary for agreement with the X-ray analysis ( il4.3.2|l . The blue 
dot-dashed line with shaded region s hows the co s molog ical simu- 
lations of metal-poor GCs from Priet o Ik Gnediiil 112003 . see i|5.1l l. 
scaled to the virial radius of NGC 1407. 

Now to find the best-fit fiducial model, we combine 
a fit of the model to the v^MsiRp) data (excluding 
the innermost point which we regard as questionable — 
see and ^3.3^ with a fit to Eq. [5] including its 
scatter. This fitting method is equivalent to a Bayesian 
analysis where the model priors are given by Eq. [5] We 
find a virial mass Mvir — 1-4 x IO^^Mq, with a corre- 
sponding virial radius of r^ir =1.0-2.0 Mpc — which is 
somewhat larger than the optical size of the group of 
0.4-0.6 Mpc (B+06b). As seen in Fig. [H (with best 
fit reported in Table [3] as model GR) , this "standard" 
ACDM model is roughly consistent with the GCS data 
(including the central wrms point) but it does not match 
particularly well the constant vyims profile outside of the 
center (—0.3 ^ 7p ^ + 0.2). The model predicts a fairly 
steep decline of ?;rms(-Rp) (slope 7p ~ +0.3), despite a 
rising Vc(r). This behavior is due partially to the strong 
radial anisotropy, and partially to the very low Sersic in- 
dex of the GCS density profile, whose flat core and sharp 
outer edge mimic the effects of radial anisotropy. 

4.1.4. Model alternatives 

Suspecting that at least one of the modeling assump- 
tions is awry, we next try excluding any a priori con- 
straint on the concentration. Without this constraint, 
there is a severe degeneracy between the virial quantities 
Mvir and Cvir, since even with an NEW profile assumed, 
we cannot make reliable extrapolations from data ex- 
tending to only ^ 5% of r^ii- As seen in the left panel 
of Fig. [Ml the DM halo could have a high mass and a 
low concentration, or alternatively a low mass and a high 
concentration. The right panel illustrates some physical 



quantities more directly constrained by the data: the 
circular velocity inside 10', and the characteristic accel- 
eration of the halo = 0.193v^/rs (which corresponds 
to a constraint on the projected mass density Psrs, or 
equivalently to the gravitational lensing convergence k); 
these quantities are more tig htly constrained t han other 
halo parameterizations (e.g. JAlam et al.|l2002f) . The GC 
data naturally prefer very low-concentration, high-mass 
halos in order to try to reproduce the flat wrms(^p) 
profile despite the low Sersic index and the presence of 
radial anisotropy. The best-fit solution is implausible 
(Mvir ~ IO^^Mq, Cvir ~ 10"^), and the fit is still not 
very good (xgc — ^•'^ ^ d.o.f.). 

Another possibility is that our adopted f]{r) profile 
is too extreme, and so we return to the standard Cvir- 
Mvir correlation, with a simple f3 — isotropic assump- 
tion. There is little theoretical justification for isotropy 
on scales of ~ 50 kpc, but there are empirical indica- 
tions of isotropic GC systems in other massive group- 
central ellipticals (e.g.. iRomanowskv fc Kochanekll2001l . 
herea fter RKOl; C+01; C+03; ISamurovic k Cirkovid 
120081 iJohnson et al.ll2009l . hereafter J+09; see summary 
in iHwang et al.l l2008l liereafter H-h08). With isotropy, 
we find a better fit to the data, with naturally a lower 
Mvir preferred since mass is no longer "hidden" by radial 
anisotropy (see model GI in Tableland in Figs. fT2l and 
[T4|) . The resulting virial radius of 0.8-1.7 Mpc is closer 
to the optical group radius. 

Given that the isotropic model provides an improved 
fit, we use this simplifying fact as a starting point to 
derive an estimated mass profile directly from the data, 
using no theoretical priors such as Eq. |31 We wish to de- 
project the wrms(^^p) profile and use it in the isotropic 
Jeans equation to find the Vc{r) profile. We start with 
the standard Abel deprojection equation for velocity dis- 
persion fe.g.. iBinnev fc Tremaind Il987l . eq. 4-58b) and 
differentiate it, integrating by parts and using Leibniz's 
Rule to find: 



r d{val) 



V 
■KV 



dr 



Rl 



dR^ 



1 d{Y.v 



2 ^ 
RMS 7 



i?D di?r, 



(7) 



where S(i?p) and j/(r) are the projected and deprojected 
density profiles of the GC system, as described in § 14.1.11 
Now parameterizing the observed fRMs(^p) by Eq. [2] 
with (wRMS.o = 225 km s"\ i?o = 4.24', 7p = -0.07; 
see Fig. 1121) . we find a steeply rising mass profile, shown 
as model II in Fig. [151 This model is similar to but 
steeper than the isotropic fiducial ACDM-based model, 
which does not reproduce the observed increasing be- 
havior of wrms(^p)- Also shown in Fig. [151 is the mass 
profile if the central dispersion points were not excluded 
(wRMS,o = 234 km s"\ i?o = 3.78', 7p = 0.24): the cen- 
tral mass would be significantly higher, with a density 
cusp inside ~ 5 kpc. 

Given that the fiducial anisotropy profile does not seem 
to be correct, and isotropy is an arbitrary assumption, 
we next investigate what the data themselves can tell 
us about the anisotropy. This transpires in two ways: 
by considering the metal-poor and metal-rich GC popu- 
lations as independent dynamical subsystems that must 
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Fig. 14. — Mass profile solutions for NGC 1407. Thick and thin solid contours show the 1 a and 2 a confidence levels, respectively (as 
found by Ax^ tests); none of the contoured solutions includes a concentration prior. The orange contours toward the bottom of the panels 
correspond to the tangentially anisotropic GCS solutions, the purple contours toward the top to the X-ray solutions (model XR), and the 
blue near-vertical contours to the group galaxies solutions. The mean relation and scatter from cosmological simulations are shown by the 
dashed curve and shaded region (see text). Additional letters label various best- fit solutions as described in the text and in Table |3] The 
left panel shows the virial mass and concentration, and the right panel shows the circular velocity inside the GC-constrained region of 10' 
(60 kpc) as well as the characteristic halo acceleration. The dot-dashed line in the right panel shows the physically-meaningful limit of the 
NFW model (i.e. Cvir = 0). The plots demonstrate that the GCS dynamics are consistent with a "normal" high-mass ACDM halo, while 
the X-ray based mass model is mildly inconsistent with ACDM and strongly inconsistent with the GC results. 



yield consistent results for the mass profile; and by mak- 
ing use of the observed kurtosis to directly constrain the 
anisotropy. 

We begin by considering the system to be scale free, 
i.e., iy(r) and (T^(r) vary with radius as simple power- 
laws with exponents —a and —7 respectively, while (3 is 
constant with radius. Then we can connect the circular 
velocity and projected dispersion by a constant relation: 

Vc{r = i?p) - kap{Rp), (8) 

where k is a com plicated function of a, /3, and 7 
(|Dekel et al.ll2005l ). We select a fixed intermediate ra- 
dius r = Rp = 4.24', where aup — 2.7, aMR — 3.4, 
7mp — —0.4, and 7mr — +0.4. Now using either the 
metal-poor or metal-rich GCs for the right-hand side of 
equation [8] must yield the same answer for the left-hand 
side, so we search for combinations of (/3mp,/3mr) that 
satisfy the equation kup/kMR — fp^MR/cp.MP- We find 
that a fully isotropic solution (/3mp = Pmr = 0) is not 
preferred, because the steeper density a^c? dispersion pro- 
files of the metal-rich GCs should depress CTp_MR relative 
to CTp^MP even more than is observed. In fact, strong 
tangential anisotropy (/? ^ — 1) is implied for both the 
metal-poor and metal-rich GCs. The degeneracies and 
uncertainties among the parameters do not allow us to 
determine more than this (if we assume /?mp — /3mr for 
simplicity, then these formally take the value of —3.6). 

We now turn to constraints from the kurtosis, which 
as higher-order LOSVD moments can provide direct in- 
formation about the orbital types. It is beyond the 
scope of this paper to explore this theme in detail, but 
we can take adva.ntage of one simplification discussed in 
iNapolitano et all (|2009l eqn BIO). If 7 = 0, /3 is a con- 



stant, and a simplified distribution function is adopted, 
then P can be directly estimated by deprojecting Kp and 
Ep — with no dynamical modeling necessary. We there- 
fore take this approach for the overall GCS of NGC 1407, 
since 7p = -0.07 ± 0.26. With Kp = -0.30 ± 0.41 (SH), 
we find that (3 = — 0.5^° |, where the uncertainty is based 
only on the observed uncertainty in Kp, and does not in- 
clude the uncertainties in Sp nor in the deviations of the 
galaxy from the simplified (7, (3) model. Thus we find an 
independent indication that the GC system is tangen- 
tially anisotropic overall. 

Given the two separate indicators of tangential 
anisotropy, we now fit a ACDM model to the GC data 
with 13 = -0.5 assumed. The best-fit model ("GT") pro- 
vides a slightly better fit than the isotropic model (see 
Fig. HH and Table [3]), and we therefore have three differ- 
ent lines of evidence for tangential anisotropy — none of 
which is statistically significant on its own, but in com- 
bination are suggestive of a reliable result. This solu- 
tion implies a slightly lower virial mass of ^ 6 x IQ^^Mq 
fFigs. [HI and [T5|) than the isotropic model. 

We can next take this best-guess model as a starting 
point to estimate the anisotropy profiles of the metal- 
poor and me tal-rich subsystems. Using the method of 
iTonrvl ()1983l . Appendix),' we take the model Vc{r) and 
the observed ap{Rp) (parameterized in our case by a 
power-law fit), then invert them to an intrinsic crr(f), 
and finally solve the Jean s equation to find /3{r) (see also 
iHobbs fc Willmorelll997l) . 

This inverse approach can be sensitive to assumptions 
about Vc and dp outside the regions constrained by the 
data, so we use it as a starting guess for simplified P(r) 
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Fig. 15.— Mass modeling results for NGC 1407 from GCS 
dynamics, expressed as circular velocity profiles. Individual models 
are labeled by letters and are described in the text and in Table |3] 
The orange shaded region around model II (direct mass inversion) 
reflects the 68% statistical uncertainties on the dispersion profile. 
The orange dashed line shows model II where the dispersion data 
at all radii have been included. 



functions, re-solving the Jeans equations and comparing 
the model cTp to the data. We adopt the simplest as- 
sumption of constant anisotropy, and find that /3mp ^ —4 
(tangentially-biased orbits) provides a good description 
for the metal-poor GCs (see Fig. [TB|) . The metal-rich 
GCs on the other hand appear to be roughly isotropic 
(/3mr ~ 0) outside ~ 2'. The high metal-rich veloc- 
ity dispersion observed at smaller radii suggests the or- 
bits become very radial in this region — although even 
/3mr = is not quite enough to reproduce the data, and 
to compensate, very tangential orbits would be needed in 
a transition region at ~ 1.5'. These peculiarities could 
be caused by an admixture of DGTOs concentrated in 
the central regions as discussed in t j3.31 and will require 
more detailed modeling in the future. Note that chang- 
ing our assumed distance and T, to higher values would 
have little impact on the inferred values for /3, Mvir, and 

Cvir ■ 

In summary, we can conclude that the GC kinematical 
information is compatible with a standard ACDM halo 
with size and mass parameters similar to a very high- 
mass galaxy group. The GCs appear to reside on some- 
what tangential orbits overall, which may be decomposed 
into very-tangential orbits for the metal-poor subpopu- 
lation, and isotropic orbits for the metal-rich GCs. The 
anisotropy results are not yet robust, and in ij4.3l we will 
explore additional constraints that may clarify the mass 
profile of the NGC 1407 system. 

For now, assuming this anisotropy result is correct, we 
report a range for the mass inside the GCS-constrained 
region of 10' (60 kpc) to be (2.7-3.2) x IO^^Mq, implying 
a mass-to-light ratio relative to the galaxy NGC 1407 of 
Ti = (29-34) Tqj, or Tb = (61-71) Tq^b- The system- 
atic uncertainty from plausible ranges on the anisotropy 
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Fig. 16.— Projected RMS velocities of NGC 1407 GCs, where 
curves show model predictions and points with error bars show the 
data. The blue dotted curve and open points show the metal-poor 
GCs, while the red solid curve and filled points show the metal-rich 
GCs. The model is the best-fit overall GC solution ("GT"), with 
/3mp = -4 and /3mr = 0. 

is ^ 20%; the uncertainties from the GCS density profile 
and the distance are at the level of ~ 10% and ~ 20%, 
respectively. The GCs with the lowest binding energies 
have apocenters of ^ 15' ^ 90 kpc, which does not nec- 
essarily mean that the GCS ends beyond this radius (cf. 
Fig. [TT|) . but rather is consistent with our model of the 
GCs on near-circular orbits. 

4.2. X-ray analysis 

We next consider a fully independent constraint on 
the mass profile of the NGC 1407 group, making use 
of X-ray emission from hot gas trapped in its poten- 
tial well. As a first approximation to estimating the 
group mass in this way, we take the iZ05'^T'-derived X- 
ray temperature of 1.02 keV (OP04), an d use t he mass- 
temperature relation of IVikhlinin et al.l (|2006( ) to esti- 
mate M500 — 3.3 X 1O^^M0 (i.e. for an overdensity of 
A = 5 00). An alternative relation from [Sanderson et al.l 
(|2006f ) that may be more valid in the group regime yields 
M500 = 1.7 X 10"Mq. With reasonable NFW extrap- 
olations to the virial radius, we thus infer Mvir ~ (2- 
6) X 10^3^0. 

A firmer result on the mass profile requires a careful 
analysis based on high-res olution X-ray data. N GC 1407 
was observed by Chandra (jWeisskopf et al.l[2002l ) on 2000 
August 16 for ~49 ksec (Obsid 791) using the ACIS in- 
strument with the S3 chip at the focus. A new level 2 
events file was produced from the level 1 events file using 
version 3.4 of the Chandra Interactive Analysis of Obser- 
vations^^ (ciAO), with CALDB version 3.3.0. Bad pixels, 
and events with ASCA grades of 1, 5, and 7, were removed 
from the analysis, and the appropriate gain map and 
time-dependent gain correction were applied. We elimi- 
nated flares by extracting lightcurves from CCDs 5 and 

^'^ http://cxc.harvard.edu/ciao/ 
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Fig. 17. — Deprojected properties of the hot gas in NGC 1407, from Chandra X-ray analysis. The points with error bars show the 
measurements in bins, where the horizontal bars show the bin boundaries in radius, and the vertical bars show the 68% uncertainties 
recovered from 100 Monte Carlo realizations of the procedure. Left: Temperature profile, computed in coarse bins, with a curve showing 
the smoothed interpolation to the binned profile. Although the fitting was done in log-log space, the plot is linear for more direct comparison 
with the GC dispersion profile (e.g.. Fig. I12I I. Right: Gas density profile, computed in fine bins, with a curve showing the cusped beta 
model. 



7, and one from the front-illuminated CCDs (2,3,6 and 
8), and filtered these lightcurves using the Markevitch 
script lc_clean. The net remaining exposure time was 
30.4 ksec. Background subtraction was performed using 
the appropriate blank-sky data set^^, which was rescaled 
using the ratio of count rates at particle-dominated ener- 
gies. We used the CIAO tool wavdetect to identify point 
sources in the data, and spectra and response files were 
extracted from the cleaned data using the CiAO analysis 
threads. 

The X-ray data analysis techniques are described in de- 
tail in J-l-09; a summary of the essential features is given 
here. Our aim is to produce a high-resolution mass pro- 
file, which is achieved through a two-stage analysis, by 
firstly analyzing temperature on a coarse grid, and then 
computing gas density on a fine grid. We extract a series 
of spectra from coarsely spaced concentric annuli, whose 
spacings are determined by a requirement of 1500 net 
counts in each annulus. We find that this criterion is 
sufficient to robustly constrain the temperature and to 
perform a successful deprojection of the spectra. We fit 
absorbed APEC models to the spectra from each coarsely 
spaced annulus using XSPEC version 11.3.2t, ignoring en- 
ergies below 0.5 keV and above 7.0 keV. We allow for 
the presence of unresolved point sources in the emission 
by including a power-law component of fixed index = 
1.56 (H rwin et ani2003[ ) following a de Vaucouleurs B}^^ 
density law with i?off= 1-17' (Z-l-07). We use the pro- 
JCT model in xSPEC under the assumption of spherical 
symmetry to deproject the spectra, fitting a single global 
abundance value across all radial bins for the sake of de- 
projection stability. This yields three-dimensional model 
parameters at a series of characteristic coarse radii. We 

http:/ /cxc. harvard. odu/contrib/maxim/bg/ 



fit smoothing spline models to the profiles using the 
SMOOTH. SPLINE algorithm from the R PROJECT statis- 
tical package^^. This allows the interpolation of the fit- 
ted parameters at a series of characteristic fine radii, de- 
scribed below. 

It is worth noting that the temperature profile 
of NGC 1407 shows typical cool-core behavior, and 
we would expect an associated abundance gradient 
(jRasmussen fc Ponmanl(2007 '). In our current modeling, 
we fit the abundance at a constant level, which affects 
the APEC model normalization and thus the gas density, 
owing to the dominance of line emission at low tempera- 
tures. The primary effect of allowing for a gradient might 
be to decrease the inferred mass at the level of ^ 20% 
(J+09). 

The resulting temperature profile is shown in Fig. [17] 
(left), where the vertical error bars are recovered from 
100 Monte Carlo realizations of the procedure, and are 
used to weight the smoothing spline fit (by the inverse 
square of the temperature uncertainty). The most no- 
table feature is a sharp peak of ~ 1.6 keV in the binned 
profile between 2' and 3'. Although this feature appears 
significant, it may be an artifact of some instability in the 
deprojection procedure; comparable analyses by II+06b 
and by Z-l-07 found a peak of only ^ 1.3 keV. The large 
error bar means that this point is down-weighted in the 
fit, but in any case, our Monte Carlo simulations do take 
large Tx excursions into account in the reported mass 
profile uncertainties. 

We next extract a series of spectra from finely spaced 
concentric annuli, and again fit absorbed APEC models 
in each annulus. We use the projct model in xspec 
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to yield the deprojected model parameters. However, in 
this stage we only fit for the APEC model normalization, 
which directly yields the gas density. The temperature, 
abundance and Galactic column in each shell are fixed at 
the interpolated value from the smoothing spline fits in 
the coarse stage. This process allows us to measure the 
gas density on a finer radial scale. 

To determine the mass profile of NGC 1407, we use 
the smoothing spline function fitted to the coarse 3D 
temperature profile, and fit a cusped "beta model" (e.g., 
H+06b) to the fine 3D gas density profile (rcoro = 33.0", 
/3x — 0.5; see Fig. [T71 right). This allows the determina- 
tion of the circular velocity profile under the assumption 
of hydrostatic equilibrium, using the equation: 

yl(^r) = ^BTx(r) / d In ^ d In Tx \ ^ 
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where Tx{r) and ng(r) are the gas temperature and den- 
sity profiles (|Fabricant et al.|[l980[ ). 

We use a Monte Carlo procedure to produce 100 real- 
izations of the mass profile of NGC 1407, which are each 
analyzed as described above. The span of results from 
these random realizations is used to define the confidence 
regions shown in Fig. 1181 where we now adopt a logarith- 
mic radius scale in order to see the models more clearly. 
We find the Vc{r) profile to be consistent with constancy 
inside 4 kpc, and outside of this radius, it rapidly rises 
until 20 kpc — increasing the enclosed mass by a factor 
of 10-15 despite a radius increase of only a factor of 5 — 
whereafter it may level off again. It is the rapid increase 
of Tx with radius (Fig.[T7]) that produces the Vc increase; 
reinstating the highest Tx point in the modeling would 
produce an even steeper (and possibly unphysical) Vc 
profile. We will further discuss the features of the mass 
profile in the next section. 

4.3. Comparing independent results 

Having derived mass profiles from both GC and X- 
ray data, we now explore the systematic uncertainties by 
making independent modeling comparisons. In §4.3.1l we 
compare our own X-ray result for NGC 1407 to literature 
X-ray studies of this system, and in ij4.3.2l we consider 
three different mass probes in NGC 1407. In i )4.3.3l we 
survey similar multiple mass probes in other systems. 

4.3.1. Multiple X-ray studies m NGC 1407 

Several independent analyses of the same Chandra 
data for NGC 1407 have been publishe d, whose results 
we also plot in Fig. [18] for reference. iFukazawa et al.l 
([2006) performed both a non-parametric deprojection 
analysis, and a parameterized "double-beta model" fit, 
to produce their mass profiles. H-|-06b also made de- 
projected models, as well as forward fitting of ACDM- 
based galaxy -|-halo models; their results extend to larger 
radii than ours by analyzing all of the available ACIS 
chips, while we conservatively used the back-illuminated 
S3 chip only. Z-l-07 again used the Chandra data, com- 
bining them also with ROSAT data, at larger radii, using 
both projected and deprojected double-beta models. 

As seen in Fig. [TH these various analyses yield qualita- 
tively similar results for the mass profile, finding a rapid 
rise in Vc{r) followed by a probable leveling-off. But in 
quantitative detail, there are differences in the results 
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Fig. 18. — X-ray-based mass modeling results for NGC 1407. 
Our new C/iandra-based result is shown as a purple solid curve with 
shaded regions showing the 1 <j and 2 a uncertainty regions. The 
best-fit stars-|-NFW solution (model XR) is shown as a pur ple dot- 
short- dashed curve. The Chandra-hsuSed results of Fukaza wa et al.l 
II2006I) are shown, including their double-beta model fit {blue dot- 
long-dashed curve) and their deprojection analysis {blue boxes with 
error bars). The Chandra-hased results of H+06b include their 
ACDM-based model fit without adiabatic contraction {red short 
dashed curve with shaded uncertainty region) and their deprojec- 
tion analysis {red circles with error bars). The Chandra+ROSAT 
analysis of Z-l-07 is shown by the green long- dashed curve, with 
green dotted curves showing the 90% uncertainties. 



that are larger than the claimed uncertainties. The Z-l-07 
results agree well with ours outside 1', but inside that ra- 
dius, have an extra "hump" in that is not seen in any 
other analysis. The hump is produced by their double- 
beta model for the gas density, a parameterization for 
which we see no motivation from our higher-resolution 
density profile. The H-|-06b Vc{r) deprojected result has 
a similar shape and amplitude to ours, but is oddly off- 
set by a factor of 2 in radius — perhaps owing to the 
different abundance treatments. We will examine these 
differences in more detail in another paper, and here note 
that there can be substantial systematic effects in de- 
riving mass profiles from X-ray data, depending on the 
parameterizations and assumptions in the modeling. 

Since the rapid Vc{r) rise between 5 and 20 kpc seems 
to be a robust conclusion from the X-ray models, we 
next consider the implications of this feature. Similar 
"kinks" have been found in the X-ray-based mass pro- 
files of many other galaxy groups and clusters (see Z-l-07 
for a summary), suggesting a transition between a DM 
halo belonging to the central galaxy, and a larger-scale 
halo pertaining to the entire group or cluster. Such a 
galaxy-group halo transition would pose a problem to 
CDM theory, which predicts that halos should relax into 
smooth mass profiles similar to the NFW profile; any 
kink would be unstable and transient. Indeed, H-|-06b 
did not find a good fit for their ACDM-based models, 
which as seen in Fig. [TH| is probably due to an inability 
of such models to reproduce the kink. The model so- 
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lutions that come closest are forced to higher DM halo 
concentrations than expected in ACDM (H+06b Fig. 4). 

We next try to fit our own galaxy+halo model se- 
quence, using our own X-ray data analysis. In contrast 
to the forward-modeling performed by H+06b, we use a 
less direct method of fitting the model Vc{r) to our de- 
projected X-ray profile. The fit is performed at 31 points 
in radius, but since there were effectively 7 independent 
data points in the original coarse-binned data, we renor- 
malize the statistic by a factor of 7/31. Our best-fit 
solution (model XR; see Table [3]) is shown in Fig. [181 
Again, the NFW-based model has difficulty reproducing 
the kink, and requires a DM halo that is more concen- 
trated than is typical for ACDM (see Fig. [H]). 

To see if we can find a better solution, we re-examine 
our modeling assumptions. The distance to the galaxy, 
and T:,, are both key parameters that could plausi- 
bly have higher values (see ii4.1.ip . However, increas- 
ing either of them causes the X-ray fit to deteriorate 
rapidly (driven by the constraints in the central re- 
gions). Decreasing either of these parameters (e.g., to 
T*j = 1.52 Tqj) yields better fits, but also pushes the 
halo concentration to higher values. 

As a single case, the high concentration of NGC 1407 
is not a problem, but we note that most of 
the X-ray bright galaxies and gr o ups e xamined by 
H+06b and by iGastaldello et al.l ()2007n were also 
found to have high concentrations, and unconvention- 
ally low va lues for T^, (se e also the high concen- 
trations in [Sc hmidt fc Allen' W?; 'DuflFv et al.' '200^; 
iBroadhursfeTa l. 2008; Broadhurst & Barkana 2008). It 
has been shown that improper modeling of the stellar 
mass in galaxies can strongly skew X-ray based infer - 
ences about their DM halos (jMamon fc Lokasl[20n5ah . 
However, in the cases of NGC 1407 as well as the other 
galaxies studied by H+06b (see also H+08), it seems that 
the X-ray data simply do not allow for plausible val- 
ues. H-|-06b found best-fits for NGC 1407 corresponding 
to T*,B ^ 2-3 Tq_b, as compared to their stellar pop- 
ulation estimate of ^ 5-6 Tq,b and ours of T*^b ~ 4- 
7 Tq_b (see §4.1.11 using our shorter adopted distance 
would lessen this discrepancy). 

We illustrate the situation further in Fig. \T9\ using the 
deprojected X-ray results of H-|-06b for their entire sam- 
ple (revising the distance to 21 Mpc for NGC 1407). 
Presented as circular velocity profiles, several trends 
are immediately apparent. Three galaxies (NGC 720, 
NGC 4125, and NGC 6482) have flat or dechning profiles 
consistent with galaxy-scale halos (modulo some issues 
with the innermost data points). The other four galax- 
ies have increasing profiles suggesti ng group-scale halos 
(cf. also' Nagino fc Matsushitall2009D . However, these in- 
creases generally set in much too quickly (by ~ 1 Ros), as 
can be seen by comparison with the ACDM-based model 
for NGC 1407. In order to reproduce these kinks, a com- 
bination of low T, and high Cvir is necessary (though still 
not sufhcient). 

The combination of awkward mass profile kinks, and 
non-standard stellar and DM mass parameters, leads us 
to suspect a systematic problem with some X-ray mass 
inferences. The kinks found by H-|-06b seem to coincide 
with strong increases in Tx, which are examples of the 
common "cool cor e" phenornenon found in galaxy groups 
and clusters (e.g. JSun et al.l[2009l ). Independently of the 
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Fig. 19. — Circular velocity profiles of seven X-ray bright galax- 
ies from H+06b. Also shown is our best-fit GC-based solution for 
NGC 1407 {solid orange curve; model GT). The open circles mark 
1 Rg.ff (in the Ks-hand) for each galaxy. Error bars show the typi- 
cal radial binning and mass uncertainty at one measurement point. 

wider implications of the cool cores, we wish to know 
whether there is a problem with applying equation [9l in 
these regions. In the case of NGC 1407, the onset of 
the kink coincides not only with a rapid Tx increase, 
but also with an onset of increasing gas isophote ellip- 
ticity (from e = 0.0 at 0.6' to e = 0.3 at 1.7') and shght 
asymmetry. It thus seems plausible that the gas in the 
kink region is not in hydrostatic equilibrium, or that 
there is additional non-thermal pressure support inside 
the kink radius (which might produce mass underesti- 
mates in these regions). Testing these possibilities with 
independent mass tracers is exactly one of the motiva- 
tions of the present GCS kinematics study, as we will see 
in the following sections. 

4.3.2. Multiple mass probes in NGC 1407 

We next compare results from the available indepen- 
dent probes on the mass profile of the NGC 1407 group. 
These constraints include the GC and X-ray analyses 
that we presented in §4. II and §4.2) as well as the dynam- 
ics of galaxies in the group (B+06b), and are summarized 
in Figs. [Hand [201 

B+06b used a friends-of-friends algorithm to identify 
19 galaxies belonging to the group, and a virial estima- 
tor to convert the kinematics into a mass. Their galaxy 
sample notably excluded NGC 1400, the bright early- 
type apparently in close three-dimensional proximity to 
NGC 1407, but with a high pecuhar velocity. B+06b 
made an alternative mass estimate including NGC 1400 
and another 5 galaxies in the region of the group — 
yielding both a higher group mass and a higher group 
luminosity; B-|-06a provided yet another updated mass 
that is intermediate to the previous two. The results for 
both the low- and high-mass group samples are shown in 
Fig. [201 implying a mass at 0.5 Mpc of (4-10) x lO^^M©, 
or Ts = 600-1200 Tq^b (the corresponding group veloc- 
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Fig. 20. — Circular velocity profile of NGC 1407, plotted vs. log 
radius. Individual models are labeled by letters and are described 
i n the text and in Table [3] Further included is our X-ray result 
( i|4.2l l. The boxes with error bars show estimates from group galaxy 
dynamics {B+06b), using the friends-of- friends sample {filled box) 
or the entire sample of galaxies in the projected region {open box). 

ity dis persion is Cp = 320-450 km s"^). iTrentham eFall 
(|2006f ) identified 10 new dwarf galaxies belonging to the 
group, and estimated a virial mass of 6.1 x IQi^^Mq. We 
plan to revisit this calculation in the future with the ad- 
dition of even more group dwarfs to the catalog. Using 
the B+06b low-mass results for now, the permitted re- 
gions of A/vir — Cvir parameter space are shown in Fig. 1141 
The orbit anisotropy of the group galaxies is an 
additional source of systematic uncertainty as with 
the GC dynamics; the virial mass estimator used by 
B-f06b implicitly assumes isotropy. Somewhat conflict- 
ing conclusions have been made about the accuracy of 
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B-|-06b found a near-zero skewness for the galaxy veloci- 
ties as well as a crossing time of only 0.02 support- 
ing the assumption of the group as virialized, and some 
indication of negative kurtosis, suggesting that isotropic 
estimators might overestimate the mass. In any case, the 
systematic mass uncertainty is probably no larger than 
the ±40% statistical uncertainty. 

Now comparing the results from the three mass probes, 
we see that they each independently imply a virial mass 
for the NGC 1407 group within the range Mvir ~ (4- 
10)xlO^^Mo (Fig. [T3]). This is an encouraging broad 
consistency, but we arc further interested in the accu- 
racy of the spatially-resolved mass profiles. The group 
galaxies and GC dynamics results are consistent (assum- 
ing fiducial DM models and any of the GC anisotropy 
profiles discussed in i i4.1.3l and §4.1.4p , providing no ev- 
idence so far of a major problem with either approach. 
But snags appear when considering the X-ray mass pro- 
file: while the virial mass from the group galaxies and 



the X-rays is compatible (Fig. [T4|) . the mass profile at 
smaller radii conflicts both with the theoretical expec- 
tations for a ACDM halo (see also i J4.3.ip and with the 
GC-derived profile. The best-fitting CDM halo param- 
eters from the GC and X-ray data are inconsistent at 
the 2 a level. This discrepancy also exists independently 
of the mass model assumed, as seen in the overlap re- 
gion (~ 2.5'-5') of Fig. [20l where the unparameterized 
mass profile inferred from the GCs is ^ 70% lower than 
the X-ray result. This difference is approximately halved 
when considering instead the deprojected X-ray profile of 
II+06b, but the problem persists. 

Our GC results are not yet robust, with the orbital 
anisotropy as the key remaining source of potential error. 
To see if this systematic uncertainty could account for the 
GCS/X-ray discrepancy, we construct a toy model where 
the X-ray-based mass profile is assumed to be correct, 
and investigate the GCS dynamics in the corresponding 
potential. To represent the X-ray result, we construct a 
Vc{r) profile with an arbitrary functional form: 



Voir) = uo + 



Vir 



(10) 



Model XT/XG has parameters 



1 + (rc/r)'; ■ 

vq — 312 km s~^, 
vi = 404 km s-\ Tc = 1.94', ( = 3.1, 6 ^ -0.37) and 
is shown in Fig. [20l the range of plausible values for the 
asymptotic profile exponent S is between —0.5 and -1-0.08. 

Assuming f3 = —0.5 for the GCS as derived in i i4.1.4l 
we carry out a Jeans analysis in this mass potential, and 
illustrate the resulting prediction for (7p(i?p) in Fig. [T^ 
(model XT). This model is a poor fit to the data, mainly 
because the relatively high X-ray-derived Vc at ~ 3' over- 
predicts CTp for the GCs around this intermediate radius. 
Setting aside our estimate for the anisotropy, we next 
investigate what profile P(r) would be required to repro- 
duce th£_GG_obserrati^^ the X-ray mass model 
(cf. iMathews fc Brighentill2003f ). Some qualitative ex- 
perimentation reveals that in order to suppress Cp at 
intermediate radii while keeping (jp high at larger radii, 
the GC orbits must be very radial in the center, and then 
rapidly become tangential at large radii. More quanti- 
tatively, we use the Tonry inversion method (see t j4.1.4p 
to derive /3(r) for the GCS, given the observed ap{Rp) 
(outside the central regions) and the assumed Vc{r). 

The result of this mass-anisotropy inversion is shown 
in Fig. 1131 the GC orbits are near-isotropic in a narrow 
range around r ~ 8', becoming very tangential at larger 
radii, and very radial at smaller radii (with an unphysi- 
cal /? > 1 for r ^ 3'). Interestingly, while the dispersion 
data inside 1.75' were not used in this analysis, the model 
does fit them a posteriori (see model XG in Fig. 1121 which 
uses a parameterized approximation to the numerically 
inverted anisotropy profile). This agreement lends cre- 
dence to the existence of some objects (possibly DGTOs; 
t i3.3p on very radial orbits near the galaxy center, but it 
is not clear that this scenario would solve the GC/X-ray 
inconsistency problem (q.v. dashed curve in Fig. I15p . 

These anisotropy conclusions are not substantially af- 
fected by the extrapolation of the X-ray mass profile (pa- 
rameter S in Eq. llOp . nor by the uncertainties in the 
GC density distribution i^(r); they are also completely 
independent of the distance assumed. Adopting the al- 
ternative X-ray based mass profile from II+06b (which 
notably accounts for a metallicity gradient) does not sub- 



Mapping the Dark Side with DEIMOS 



21 



stantiaUy change the /3{r) result. Reconcihng the GC 
and X-ray data therefore seems to require an ad hoc GC 
anisotropy profile with extreme radial variations, con- 
flicting with the anisotropy inferred from the GC data 
themselves. As we will discuss in §5.11 this profile is also 
highly implausible on theoretical grounds. 

We therefore have apparently found a fundamental in- 
compatibility between the GC and X-ray results. How- 
ever, it is premature at this juncture to pinpoint the 
culprit, since there are puzzling quirks in both analyses 
(the X-ray kink and the central GC dispersion peak). 
More exhaustive study of this system will be needed in 
the future, including stronger empirical constraints on 
the GC anisotropy (see t j4.3.3p . and a dynamical analy- 
sis of the field star kinematics — which should provide a 
crucial cross-check on the central X-ray mass profile (cf. 
the study of M60 by H-l-08, where the mass inside 1 i?eff 
derived from stellar dynamics may be significantly higher 
than the X-ray-based mass) . 

For now, we wish to estimate the best mass model that 
jointly fits all of the constraints, surmising in the spirit 
of the central-limit theorem that the various systematic 
problems will roughly cancel out. The inconsistencies 
between the different constraints do suggest that we not 
naively use the usual method of simply adding up all of 
their values, but rather attempt to account for the 
systematic uncertainties through a weighting scheme. In 
particular, we do not wish the tighter constraints from 
the X-ray analysis to dominate the final solution. 

Combining weighted data set s is a ii on-trivial exer- 
cise (see e.g.. Press 1997; Hob son et al.l 12002) . and for 
now we adopt the following ad hoc approach. We renor- 
malize the corresponding to the X-ray and GC con- 
straints by their best-fit values when each of them is fit- 
ted independently, including the ACDM concentration 
prior. The new functions, as well as the unmodified 
group-galaxies x^ (corresponding to the total mass in- 
side 500 kpc; see Fig. EH)) , are summed along with the 
concentration prior to produce a new total x^- 

The resulting "consensus" solution is a normal ACDM 
halo with Mvir 7 x lO^^M© (see model C in Table [3] 
and Figs. fT4l and [20)1 . With this solution, DM becomes 
dominant over the baryonic mass outside of r ^ 1.3' 
(r-^ 1.4 i?cff ~ 8 kpc). The benchmark mass-to-light ra- 
tio gradient parameter introduced by iNap olitano et all 
(11)05) is ViT - 2.0 for NGC 1407, which is by far the 
highest value reported for any early-type galaxy, demon- 
strating the extreme dominance of DM in this system. 
Note that the DM scale radius is much larger than 
the extent of the GC kinematics data (120 kpc vs 60 
kpc) , so it is really the group galaxy kinematics that de- 
termines Mvir, with Cvir then pinned down by the GC 
and X-ray constraints^^. In the near future, it should be 
possible to push the GC kinematics measurements out 
to larger radii, and gain a stronger constraint on Mvir 
without reference to the group galaxies. 

One alternative estimate for the mass profile would 
be derived from the GCs and group galaxies only, dis- 
regarding the X-ray data and the ACDM concentra- 
tion prior. The resulting model GG has a mass of 

The crossing time for the outer GCs is only ~ 0.2 Gyr, and 
for the group galaxies is ^ 0.3 Gyr, so non-equilibrium effects are 
unlikely to be the explanation for the high mass inferred. 



Afvir ~ 1.1 X 10 M0 and a slightly low concentration 
(Fig. [HI). 

4.3.3. Multiple mass probes in other systems 

More detailed analyses of the current and future data, 
including direct constraints on the GCS anisotropy, 
should provide a more accurate picture of the NGC 1407 
mass profile. But given the apparent conflicts in this 
case between independent mass constraints, we pause to 
review similar comparisons that are available in other 
early-type galaxies. Such an approach has become com- 
monplace in the arena of galaxy clusters, where gravita- 
tional lensing. X-rays, and galaxy dynamics may all be 
used as cross-checks and to provide st ronger const raints 
when used jointly (e . g.. iGirardi et al . 1998; Sand et al.l 
120041 iBenatov et all 120061 : iMahdavi et al.i ,2008. ). But 
comparable studies in galaxies are in their infancy, as 
the data are difficult to acquire. 

A handful of galaxy halos have been studied using both 
X-rays and kinematics; these are mostly cases of central 
group and cluster galaxies with GC data, b ut sometimes 
also with PN kinemati cs (e.g., ICohen &z R vzhov 199^ 
Kronawitter et al. '2000"; R KOl; C+"oi: iNapolitanq et ail 



2002t C-H03: iRichtler et' all 120041; [Ixikazawa et al.ll2006t 



Bridges et all l2006t ISchuberth et al.l l2006t ISamurovid 
20061 IWu &: Tr emaine' '2006'; 'Richtler et al.' '2008!, here- 
after R-H08:IChakrabarty & Raychaudhury 2008; J+09). 
Within the uncertainties, most of these X-ray /kinematics 
inter-comparisons are broadly consistent: the discrep- 
ancy we find in NGC 1407 is one of the largest seen^^. 
However, the kinematical tracers in these studies have 
generally been assumed to be fairly isotropic — which may 
indeed turn out to be typical for the GCs in group/cluster 
central galaxies, but is not yet certain. 

Given the freedom to adjust orbital anisotropics, 
one could make almost any two data sets agree (e.g., 
§4.3.2p . so to conclusively compare mass results, one 
must derive independent constraints on anisotropy. This 
is possible by measuring higher-order shapes of the 
LOSVDs, which can be used to constrain the types 
of orbits (cf. our crude kurtosis analysis in i i4.1.4p . 
Such mass/anisotropy degeneracy breaking techniques 
are well-established for integrated-light stellar kinemat- 
ics, bu t can als o be applied to discrete velocities (e.g. , 
RKOl: iRomano wskv et al. 2003: IWu fc Tremaind l2006l : 
iDe Lorenzi eral.,,2008^ .Chaname et al.ll2008l ). Although 
it is commonly supposed that a few hundred veloci- 
ties are not sufficient to break this degeneracy, in fact 
additional constraints may be brought to bear on the 
problem. One is to make use of stellar kinematics in- 
cluding their LOSVDs to pin down the mass profile 
out to ~ 1-2 i?cff (cf. RKOl). Another is to assume 
that the underlying spatial profile of the tracer popu- 
lation is fairly well known, e.g., the total GC surface 
density profile represents the sa me population as the 
kinematically-measured G Cs (cf. iMerritt fc Sahal[l993l : 
iMerritt fc Tremblavl[l99l . 

So far, the only galaxy to have its DM content studied 
using both high-quality X-ray and optical data (including 
LOSVD constraints) is M87. RKOl used nonparametric 
"orbit modeling" of the dynamics of the GCs and field 



^^ IMathews fc Brighentl l|2003l) found a discrepancy between X- 
ray and dynamics-based mass profiles in the central regions of M49. 
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Fig. 21. — Mass modeling results for M87. The orange solid 
curve shows the best fitted ACDM-based orbit model to the stellar 
and GCS kinematics; the dot-dashed curves show other acceptable 
solutions (RKOl). The filled circle shows the outer extent in ra- 
dius of the kinematics data fitted. The orange dashed curve shows 
the result of anisotropic distribution-function based modeling of 
the GCS kinematics by Wu & Tremainc (2006), with the extent of 
the curve showing the radial range of the data fitted. The pur- 
ple shaded region show s the XMM-Newton-ha.sed mass model of 
[Matsushita et all (1200211 . The blue jagged solid line with shaded 
region shows the Chandra-hased model of C-(-08 (mean molecu- 
lar weight fi = 0.61 assumed), with 4x radial binning applied to 
smooth out the highest- frequency variations; the uncertainty region 
incorporates the differences from single- and multi-temperature 
modeling. The black dotted curve shows a model with no DM. 
The figure thus shows that the detailed optical kinematics con- 
straints are broadly consistent with the X-ray mass results, but 
without artificial smoothing, the X-ray data imply an unphysical 
mass profile. 

stars in this galaxy, deriving a mass profile that may be 
compar ed wi th t h e sub sequent XMM-Newton result of 
[Matsushita et^all (|2002f ). As seen in Fig. 1211 the optical 
and X-ray results in the region of overlap are nicely con- 
sistent, except at small radii (<2 kpc), where the typical 
problem of an X-ray core mass underestimate is apparent 
( e.g.,fer iRhcnti fc Mathews 1997; H+08). 

IChurazov et al.l (j2008i hereafter C-l-08) have recently 
used new Chandra data to examine the mass profile of 
M87 in more detail, and again found that the optical and 
X-ray results on the mass profile are broadly consistent, 
but that the X-ray profile contains rapid radial variations 
that are not seen (by construction) in the smooth opti- 
cal models. These "wiggles" appear at a nuisance level 
when examining the gravitational potential, but when 
differentiated to derive a profile of Vc{r) or M(r), the 
wiggles are amplified catastrophically, even implying at 
one locatio n 10 kpc) a negative enclosed mass (see 
Fig.[2T]and rGhizzardi et aLlfeOOd ). The wiggles were not 
seen in the earlier XMM-Newton study because of its use 
of smooth analytical models to describe the gas density 
and temperature profiles. This standard practice avoids 
amplifying noise in \ow-S/N data, but in the case of M87, 
it appears to have suppressed real features in the data 
(see also J+09). 



C-l-08 showed that the apparent M87 mass wig- 
gles could be readily explained as a consequence of 
shock waves propagating through the gaseous interstellar 
medium. They claimed that one can ignore these local- 
ized effects and still safely consider global changes in the 
mass potential, taking also into account additional (ap- 
parently minor) corrections to the X-ray based mass such 
as non-thermal pressure support (see also'Pizzolato et alj 
2003; Ciotti & Pellegrini 2004; Pellegrini & Ciotti 200a 
Rasia et al.l l2006t Nagai ct_aL| 120071 : iDiehl fc Statld 
20071: lO'Sullivan et al. 2007). However, in practice, it 
is difficult to know a priori how to smooth out the X-ray 
wiggles, and thus the conclusions derived about DM halo 
properties in detail can be strongly affected. NGC 1407 
is a case in point, as its X-ray kink feature could be 
attributed to a shock, but because we do not have more 
extended X-ray data to provide the asymptotic mass pro- 
file, the feature forces the ACDM models toward high 
concentrations and low values for T*. However, there is 
less a priori reason to suspect AGN-generated shocks in 
NGC 1407, as its X-ray image shows no major distur- 
bances, and its radio activity is '-^ 1000 times less than 
that of M87 (iCondon et al.lflQOl . 

M87 has clearly disturbed X-ray isophotes, while there 
may also exist a class of fairly undisturbed systems where 
the X-ray results will be very reliable. JDiehl & Statle^ 
(2008b) have introduced a parameter 77 to describe asym- 
metries in gaseous X-ray halos. The correlations of rj 
with central radio and X-ray luminosity, and with outer 
temperature gradient, suggest that both the ambient 
medium and the central AGN cause disturbances in the 
X-ray gas — even in cases such as N1407 where the AGN 
activity appears fairly weak. 

Values of 77 were provided for most of the galax- 
ies with X-ray /kinematics mass comparisons (see ref- 
erences at the beginning of tj4.3.3p . The galaxies are 
M60, M49, NGC 1399, NGC 4636, and NGC 1407, 
with r] = 0.05, 0.08, < 0.09,0.09, and 0.17, respectively. 
NGC 1407 has a relatively high value of 77 (despite having 
less obvious disturbances than e.g. NGC 4636) and man- 
ifests a mass discrepancy, but there is otherwise no obvi- 
ous correlation between 77 and the reliability of the X-ray 
results. So far, none of the other galaxies clearly shows 
a large mass discrepancy, but all of them have at least a 
hint of a non-equilibrium kink in the mass profile — even 
in the case of M60, which has one o f the smoothest gas 
halos in the IDiehl fc Statled ()2008bD sample. 

On a related note, 4 of the 7 galaxies studied in X- 
rays by H-|-06b have deprojected mass profiles that are 
inconsistent with the forward-fitted smooth NFW-based 
models (see i i4.3.ip . These are again cases where the X- 
ray approach may be problematic, but the rj values for 
them and the other 4 galaxies are not clearly different. 
C-l-08 did note that 77 may not be a reliable diagnostic of 
hydrostatic equilibrium, and we speculate that a metric 
based on radial (rather than azimuthal) fluctuations in 
the X-ray photometry and/or temperature may be more 
useful for this purpose. 

The generic reliability of X-ray mass constraints could 
be further tested by applying robust dynamical models 
(as in the case of M87) to the stellar and GC data in 
NGC 1407 and in the other galaxies mentioned (while 
also extending the models to consider non-spherical ef- 
fects). These are demanding exercises that we defer to 
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TABLE 4 

KiNEMATICAL PROPERTIES OF NEARBY GROUP/CLUSTER-CENTRAL GALAXIES. 



NGC Type Env. A^qc D Mb ctq v^ot/up frot/cp tiRMS fRMS Rcf. 
[Mpc] [km s-^] (MPGCs) (MRGCs) (MPGCs) (MRGCs) 



1399 


El 


CCG 


228 


18.5 


-21.1 


344 ±6 


IQ+O-ii 


n+0-10 




258ii7 


D+04, R+04 


1407 


EO 


GCG 


172 


20.9 


-21.0 


273 ±6 


U-18_o.l8 




234121 


247li7 


this paper 


3379 (M105) 


El 


GCG 


38 


9.8 


-19.9 


207 ±2 










B+06c 


4472 (M49) 


E2 


GCG 


127 


15.1 


-21.7 


291 ±8 




26+0-i"3 


0( 0_2g 




C+03 


4486 (M87) 


E3 


CCG 


144 


14.9 


-21.5 


332 ±5 




66+0-iS 

U.DD_(j 23 


414li 


383+^^ 


H+01 


4636 


E2 


GCG 


88 


13.6 


-19.8 


203 ±3 




'^■"'-0.57 


209^25 


180+26 


S+06 


4649 (M60) 


E2 


GCG 


59 


15.6 


-21.3 


336 ±4 


U.48_o 28 


'-'■^'^-0.53 




178l« 


L+08 


5128 


E2/S0 


GCG 


170 


3.9 


-20.7 


120 ±7 


54+0-24 


40+0-10 


75± 11 


15911^ 


W+07 



Note. — The third column denotes cluster-central galaxy or group-central galaxy. The fourth column is the number of GCs with 
measured velocities (after discarding outliers). The seventh column shows the galaxy central velocity dispersion from HyperCat". The 
eighth through eleventh columns show parameters for the metal-poor and metal-rich GCs (MPGCs, MRGCs), as reca lculated for this paper. 
The da ta references ar e H+0 1: |Hanes et al. (2 00j); C+03: Cote et al. (2003); D+04: Dirsch et al. ( 200 |); R+04: IRichtler et al.l || 2004 ): 
B+06c: IBergond et"an l l200l) : S+06: .Schuberth et al.l (|200"Q1 : W+07: iWoodlev et al., (,2007. '): L+08: iLee etaLl (I2008bll 



''http://leda.univ-lyonl.fr// dPaturel et al.ll2003r) . 

future papers. One might consider the X-ray analyses 
more suspect owing to the highly uncertain physics in- 
volved, but even in the relatively well-posed case of grav- 
itational dynamics, there are questio ns about the relia - 
bility of the modeling results (e.g., iDekel et"aI1 12005D . 
The latter issue can be addressed by cross-checking in- 
dependent kinematical constraints, e.g., from PNe, and 
from the metal-poor and metal-rich GC subsystems sep- 
arately. We hope in the future to pursue these avenues 
in NGC 1407 as well as in other galaxies; studies with 
GCs and PNe will also be useful for testing the presumed 
association between metal-rich GCs and field stars. 

5. THE GROUP CONTEXT 

Having determined some properties of the mass dis- 
tribution and GC kinematics in the NGC 1407 group, 
we now attempt to place this system in the context of 
what is known and expected for galaxy groups and clus- 
ters. Are the high T value and rising Vc{r) profile really 
unusual for groups of comparable luminosity? We may 
alternatively consider the system as a dim cluster, since 
its mass approaches the usual (arbitrary) group/cluster 
boundary of IQ^'^Mq: is it rare for clusters to have been 
so inefficient in forming stars? How do the GCS kine- 
matics compare to those of other systems? 

To address these questions, we compile all of the 
published data available for GCS kinematics in group- 
and cluster-central galaxies (generically abbreviated as 
GCGs) with reasonably large and radially-extended GC 
velocity catalogs (see Table 2]). For the sake of consis- 
tency, we re-derive the kinematical properties of these 
data sets using the same techniques that we have applied 
to NGC 1407. The redshifts of all the galaxies besides 
NGC 1407 are low enough to cause potential contami- 
nation from Galactic stars, which can easily distort the 
kinematical inferences by their accidental inclusion, or 
by over-compensating and excluding some low-velocity 
GCs. To be conservative, we have rejected all objects 
with velocities below Vsys of the host galaxy 

1^ A similar exercise for NGC 1407 yields rotation and kurtosis 
parameters that are formally inconsistent with the full data-set 
results, which is a warning that the uncertainties in the kinematical 



We discuss the GCS kinematics and dynamics in ij5.ll 
the mass profiles in H5.21 the baryon content in H5.3[ and 
some implications in i |5.4l 

5.1. GCS kinematics and dynamics 

The kinematics of any galaxy should provide clues 
about its evolutionary history. However, GCGs have 
fairly similar optical properties to their mo re isola ted el- 
liptic al galaxy counterpar t s (e.g., Ivon d"er Linde n et al.l 
I2OO7I: IBildfell et al] [20081 ILoulpser et aP 120081) . which 
might refiect an early galaxy collapse history decoupled 
from the more recently assembled surrounding group. To 
probe the assemblies of both the GCGs and their groups, 
the large radial range of the GCs could be of service. The 
properties we will consider include the rotation ( i|5.1.ip . 
velocity dispersion ( iJ5.1.2p . and orbits (' t|5.1.3p of GCs. 

5.1.1. Rotation 

We begin by examining GCS rotation in the galaxy 
sample, considering the metal-poor and the metal-rich 
rotational dominance parameter (frot/cp) and the kine- 
matic misalignment angle (^o.kin ~ ^o.phot, where the 
±180° degeneracy in 6*0, phot is broken by reference to the 
direction of stellar rotation from the literature). These 
parameters are in general very dependent on the sam- 
pling of the GC velocity data by color, magnitude, and 
position — as illustrated in Table [5] and Fig. [21] by the re- 
sults for different radial bins of NGC 1407 GCs. However, 
momentarily throwing caution to the wind, we plot all 
of the available data from Table [Ij (except for NGC 3379 
which does not have enough GC velocities for this pur- 
pose). 

These systems are generally consistent with an over- 
all mild rotation of v^ot/crp 0.2-0.6 seen equally in 
both the metal-poor and metal-rich GCs (see also dis- 
cussions m iBeaslev et alll2006t IWoodlev et alll2007l : and 
H+08). The kinematic position angles of the metal- 
poor and metal-rich GCs are in most cases consistent 
with each other, but with no preferred alignment rel- 
ative to the galaxy isophotes (i.e. both major- and 

parameters are probably underestimated in general. 



24 



Romanowsky, Strader, Spitler, Johnson, Brodie, Forbes, & Ponman 



0.8 - 



u 

1 0.6 



h°0.4 - 



0.2 - 



TT — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



J 50 

BO 

<U 

-o 



x: 
o 



ffllOO 



I 



50 



-1 1 1 1 1 1 1 1 1 1 1 — -p 1 1 1 1 1— 

1^ 1 



I 



_l I L_ 



N1407-OUI. 



_l I U L 



0.2 0.4 0.6 0.6 

(^rotAp) (metal-poor) 



50 



100 



150 



I'^o.kin-'^o.phoJ (metal-poor) [deg] 



Fig. 22. — GCS rotation parameters for group/clustcr-central galaxies, with the metal-poor and metal-rich subsystems compared. In the 
case of NGC 1407, separate points are plotted for the overall GCS, for intermediate radii (1.77'-4.52'), and for outer radii (4.55'-ll.l'). 
The dotted diagonal lines mark equal values for the subsystems. Left: Rotation amplitude relative to velocity dispersion. The X marks 
the results of the merger simulations of B-l-05. Right: Kinematic misalignment. The angles (?o,kin and So, phot correspond to the angular 
momentum vector and to the photometric minor axis, respectively. Classic rotation around the minor axis is at 0° , rotation along it at 90° 
(dashed lines), and counter-rotation at 180°. 

increases with radius. 

B-l-05 also predicted that the GCs would have kine- 
matic misalignments relative to the galaxy isophotes, 
with similar rotation axes for the metal-rich and metal- 
poor subsystems — again in broad agreement with the ob- 
servations. It is not clear whether kinematic twists with 
radius, or kinematic misalignments between the stars and 
GCs, are expected from these merger models. Intrigu- 
ingly, B-l-05 also found that the rotation of a GC system 
aligns with the rotation of its host DM halo, potentially 
providing an indirect probe of an otherwise inaccessible 
property. 

Detailed kinematics predictions do not exist for other 
GCS formation scenarios, such as accretion, quasi- 
monolithic multi-mode galaxy collapse, and metal-rich 
GC formation by major mergers — but with all of these 
scenarios, one would probably expect a difference be- 
tween the rotation of the metal-poor and metal-rich 
GCs that is not seen in the data (see discussion in 
H-l-08). This observation suggests that the GCs were 
in general formed prior to the assembly of their host 
GCG galaxies by mergers — which affected the metal- 
rich and metal-poor GCs equally, replacing any pre- 
existing orderly rotation with merger-induced hot rota- 
tion. A related implication is that the last major merger 
episodes in the se galaxies were "d ry" or "damp" (rela- 
tively gas-DOo n iForbes et al.ll2007|), as also suggested by 
other studies dP e Lucia fc Blaizotf '20071: iLiu et al.l[200a 
iMcIntosh et~aL 2008: Tran et al. 200|). 

Such a scenario of early GC formation was considered 
as part of a recent foray in to the m odeling of GCGs in 
a cosmological context (Bek ki et al.lf2008j ). Their predic- 
tion of relatively low GCS rotation overall (frot/cp ^ 0.3) 
appears to conflict with the observations, but the simula- 
tions did not directly include non-gravitational processes, 



minor-axis rotation are found, as well as orientations in- 
between). Such kinematic misalignments are unsurpris- 
ing, since the photometric axis of a t riaxial galaxy can 
be readily twisted in projection (e.g., iFranx et al.lll991l : 
iRomanowskv fc Kochanekl ll998^. 

The GCs also appear in some cases to be misaligned 
relative to the (central) stellar kinematics, even to the 
extent of counter-rotation, which would imply an intrin- 
sic kinematic misalignment. Furthermore, there are ap- 
parent kinematical twists with radius (see ^3.2]i — as in 
the case of M87, where such a twist may coincide with 
the onset of the cD envelope and reflect a transition 
to the group environment (see C-l-01). Note however 
that a stellar/GC misalignment, or a kinematical twist 
within a population, does not necessarily imply radi- 
cally decoupled sub-components, since subtle shifts in 
the orbit populations with radius c an dramatically al- 
ter the net angular mome ntum (e.g.. iJesseit et al.]l2007l : 
Ivan den Bosch et al.|[200l . 

Deciphering the implications for the intrinsic structure 
and kinematics of these galaxies and their GC systems 
will require improved data and in-depth modeling. But 
our current overall impression is that the metal-poor and 
metal-rich GC subsystems have similar rotational prop- 
erties, which might be distinct from the rotation of the 
field stars. 

Considering next the connections to galaxy formation 
theory, we plot in Fig.[22](Ze/i() a theoretical prediction for 
the rotation amplitude from B+05, as the median value 
and spread inside 6 i?off of simulations of dissipationless 
disk galaxy pair major mergers (their multiple merger 
results appear to be fairly similar). Broadly speaking, 
their predictions (t'rot/cp ~ 0.4 for both metallicity sub- 
systems) agree with the data. We will leave it as a fu- 
ture exercise to check their prediction that the rotation 
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which could alter the rotation predictions. Also, their re- 
sult reflects all of the GCs within the virial radius, while 
the observations sample the GCs within only ~ 0.1 Tvir- 

It will be interesting to see if a different rotation pat- 
tern emerges from observations of non-central ellipticals 
(e.g., B. Kumar et al., in prep.), and of the general 
class of fainter, d isky, rotation-dominat ed ellipticals and 
lenticulars (e.g., lEmsellem et al.l [2007( ). Such systems 
may have undergone more sporadic merger histories in- 
volving different mass ratios and more recent GC forma- 
tion, and may have experienced more outer halo tidal 
stripping than GCGs have. Measurements of PN kine- 
matics in all of these systems would also allow testing of 
the B-l-05 prediction that the GCs and the field stellar 
halos should have kinematical similarities (br oadly con- 
firmed in the cases of NGC 5128 and M31 ; iPeng et al.l 
[200l iWoodlev et al.l[2007l: iLee et al.l[2008ah . 

To be robust, the theoretical predictions must also be 
improved to include baryonic effects in the mergers, and 
a full cosmological context. 

5.1.2. Velocity dispersions 

We next compare the velocity dispersions in the GCG 
sample. In Fig. [221 we have applied our smoothing 
technique to the data sets listed in Table S] to produce 
GCS dispersion profiles. We note first of all that our 
NGC 1407 data out to 60 kpc comprise the most radi- 
ally extended set of GC velocities so far published for an 
early-type galaxy (a forthcoming data set for NGC 1399 
will reach 80 kpc, as previewed in R-l-08). Next we see 
that in GCGs, a fairly constant profile is so far a generic 
property (see also the discussion in H-l-08). Some depar- 
tures from constancy are also evident, with hints in sev- 
eral galaxies of a central dispersion drop to '^15-20 kpc, 
driven by the metal-rich GCs (with NGC 1407 the most 
extreme case). M87 also stands out with its rapidly as- 
cending dispersion. Although there is some indication 
in the NGC 1407 data of the onset of a large halo as 
in M87, the mass scales involved are considerably differ- 
ent (we will examine their mass profiles more directly in 



In the context of the B+05 simulations, the fairly con- 
stant GC dispersion profiles in all of these systems would 
support their formation in multiple mergers, since pair 
mergers were found to have systematically more steeply 
declining dispersion profiles. This result in the simula- 
tions is presumably caused not by mass profile differences 
but by some combination of variations in the GCS spatial 
densities and orbital anisotropics. However, none of the 
B-l-05 merger remnants shows a dispersion profile quite 
as flat as seen in the data — a discrepancy that could very 
well be due to more massive, extended DM halos around 
these GCGs. Interpreting the patterns seen in Fig. [53] 
requires dynamical models to include the effects of mass, 
anisotropy, and spatial density (e.g., H5.2\i . 

We next make a basic comparison between the global 
velocity dispersion amplitudes of the metal-poor and 
metal-rich GCs in these systems, and notice some inter- 
esting trends (Fig.[211). The urms of the metal-poor GCs 
seems to roughly track the central velocity dispersion of 
the host galaxy (ao), while metal-rich GCs may have 
a weaker dependence on ctq. One would naively expect 
the opposite, since the metal-poor GCs extend farther 
out into the group halo and should be less coupled to 
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Fig. 24. — Overall GCS velocity dispersions vs central galaxy 
velocity dispersions. Open blue circles show metal-poor subcom- 
poments, and filled red circles show metal-rich. The dotted line 
shows a one-to-one relation. 

the central galaxy's properties. Clearer conclusions will 
require a more in-depth, homogeneous analysis of these 
and additional systems. 

5.1.3. Orbital properties 

Lastly, we consider the dynamics of the GCSs, in the 
context of their orbital anisotropics. We have made a 
tentative determination for NGC 1407 that the over- 
all GCS has somewhat tangential anisotropy, which is 
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for the metal-poor and metal-rich GC subpopulations of GCGs. 
Dashed lines mark Gaussian distributions (Kp = 0). Negative and 
positive kurtosis generally imply tangential and radial anisotropy, 
respectively. 

strongest for the metal-poor GCs, with near-isotropy for 
the metal-rich GCs fi i4.1.4p . Such studies are in their in- 
fancy, but there is an emerging consensus that the overall 
GCS anisotropics of group-central ellipticals are isotropic 
or modestly tangential (see overview in H+08). To sum- 
marize the situation, we plot the kurtoses of the metal- 
poor and metal-rich systems in Fig.[25l In all cases, both 
subpopulations are consistent with zero or negative kur- 
tosis (Kp = —0.2 ± 0.2 on average), suggesting isotropic 
or somewhat tangentially anisotropic orbits. 

This observational result conflicts with the general ex- 
pectation from theory that stars and GCs in galaxy ha- 
los should be in somewhat radial orbits owing to infall 
and merger processes (see references in i j4.1.3p . a pre- 
diction that has been observat i onally confirmed for PN e 
pe Lorenzi et al.l 120081 120091 : iNapolitano et all I2009D . 
In particular, the results for the metal-poor GCs stand in 
stark contrast to the proposal of D-f 05 and iMoore et al.l 
([2006') that these GCs are associated with DM subsets 
on radial orbits (see Fig. [T5)) . 

One caveat is that the theoretical predictions have gen- 
erally considered individual galaxies rather than GCGs, 
which could in principle have distinctly different halo and 
GCS properties — e.g., perhaps the transition radius tq in 
Eq.[6]would be much larger in such systems^^. There are 
so far no comparable data on more isolated galaxies, ex- 
cept for the Milky Way, where the metal-poor GC orbits 
are o nly mildly radial on average ()van den Bosch et al.l 
[19991) . 

The other key consideration is that the present-day 
population of GCs is probably a small surviving rem- 
nant of a primordial population that was subject to 



The simulations of GCG halo stars bv lSommer- LarsenI pOOg ) 
pre dict anisotropics that are very similar to the D-l-05 model in 
Fig. 1131 f GR.XR*) in the case of fossil groups, and slightly less radial 
in the case of non- fossils. 



various destructive effects such as ev aporation, dynami- 
cal fr i ction, and tidal shoc king (e.g.. iGnedin fc Ostriked 
Il997t iFall fc Zhai^ 1200 In principle, these effects 
could be correlated with orbit type, so that GCs on 
radial orbits were preferentially destroyed, leaving a 
more isotropic or tangential distribution for present-day 
observers Ce.g.. lAeuilar et al.l 119881 : iMurali fc Weinberg 
fl997t iBaumeardtl I1998D . However, a simulation of 
GCs in situ around an M87 analog indicated that or- 
bit anisotropies a re decreased only wit hin a few kpc of 
the galaxy center (|VesDerim et al.ll2003[ ) — while the GCS 
anisotropies are observed to be low on scales of tens of 
kpc. 

Still, there is much work remaining to be done for 
the predictions of initial GC orbits, and the connections 
of orbits with destruction. More realistic treatments of 
shocks and triaxial orbits may imply that GC disruption 
is more effective at larger radii (e. g., iDehnen et aLll2004l : 
ICapuzzo-Dolcetta fc Vicarl 120051 ). Most importantly, 
the evolution of GCs should be followed in the context 
of a full cosmological picture of galax y assembly. The 
only w ork along these lines so far is bv lPrieto fc GnedinI 
()2008f ). who found that tidal shocking is not a domi- 
nant effect in metal-poor GC evolution, and that the 
GC orbits are set by the orbits of their progenitor satel- 
lite galaxies (whether extant or disrupted), resulting in 
near-isotropy out to radii of ~ 0.15 rvir, which is not 
too different from our observational constraints (Fig.[T3]). 
However, the GCSs in these simulations are much more 
spatially extended than real systems, calling into ques- 
tion their mechanism for depos iting the GC s and thus 
their anisotropy predictions (cf. lAbadi et al.i r2006) . Fur- 
thermore, the simulations describe an isolated Milky Way 
counterpart rather than a GCG, where the violent merg- 
ing process should erase whatever initial orbits the GCs 
had anyway, leaving again the radial orbit puzzle dis- 
cussed above. 

Keeping in mind that the disruption rates of GCs may 
depend on their masses and concentrations, one wonders 
whether there is any signature dependence of GC kine- 
matics with their luminosities or sizes. In fact, there 
are only a few galaxies where this issue has been ex- 
amined: NGC 1 399, N GC 4374, and NGC 4636 (R-t-04; 
ISchuberth et al.l 120061 : B. Kumar et al., in prep.), and 
now NGC 1407 (see EH). In each case, a plot of GC 
velocities with magnitude reveals that the very brightest 
objects have noticeably platykurtotic distributions and 
correspondingly higher dispersions than the rest of the 
GCs. These features are not statistically significant in 
all cases (including NGC 1407), but the overall pattern 
is clear. One explanation could be that the population 
of bright GCs on radial orbits was depleted by disrup- 
tion effects, leaving them with only near-circular orbits. 
This does not necessarily mean that the brightest GCs 
were most sensitive to these effects, since a general re- 
duction in mass for all GCs on radial orbits would shift 
them toward lower luminosities, leaving the bright end 
of the distribution as the only subpopulation observably 
depleted of radial orbits. 

In NGC 1407, we also see peculiar behavior at small 
radii, where the metal-rich GCs have a sharply increasing 
velocity dispersion that could imply a radial orbit distri- 
bution (see H3.3\ and ii4.3.2p . This conclusion seems to 
contradict the generic theoretical expectation that ccn- 
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tral GCs will tend to more tangential orbits, as high- 
velocity radially-plunging objects would be destroyed 
near their pericenters. The presence of the DGTO can- 
didates fi )3.ip suggests another possibility: that some of 
the "normal" central GCs are really DGTOs that were 
formed from the nuclei of tidally-stripped g alaxies (e.g. , 
Bekki et al.l[200H iFellhauer fc Kroupall200a lBekkill2007l : 



Goerdt et al.l 1200^ 1 These DGTOs would now have 
radially-biased orbits, while a population of tangentially- 
biased progenitor galaxies should linger at larger radii, 
for the same reasons discussed above for normal GC dis- 
ruption. Whether the detailed anisotropy profiles from 
this theoretical scenario and the NGC 1407 observations 
are consistent, and whether the progenitor population 
can be found, remains to be seen. Other open questions 
are why there would be no signature of a residual central 
DGTO population in the GC kinematics of other GCGs, 
and why the DGTO population around NGC 1399 is ob- 
served to have a decreasing ve locity dispersion toward 
the center (jThomas et al.ll2'0Q8[ ). 

5.2. Mass profiles 

We next examine the mass profiles of the systems 
from Table [3] (we have not re-derived these profiles as 
we did for the kinematics properties). For the sake 
of homogeneity, these include only results where the 
halo masses are constrained by the GCs rather than 
by any independent constraints such as PNe and X- 
rays. X-ray studies in particular could be selection- 
biased toward system s with unusuall y massi ve, concen- 
trated DM halos (e.g..lyikhlinin et al. 2006; Buote et all 
I2007t iFedeh et al.l l2007[). since the resulting high X-ray 
luminosities of such systems are tempting observational 
targets. The GC studies so far may also have a similar, 
less direct bias, since the galaxies hosting large, readily- 
studied GC pop ulations are also likel y to have more mas- 
sive DM halos (jSpitler et al.ll2008b( l. Such a bias suits 
our current analysis of GCGs, but more general conclu- 
sions about halo masses in elliptical galaxies will require 
a careful galaxy sample selection independent of the GC 
content (e.g., B. Kumar et al., in prep.). 

All of these GC-based mass results involved near- 
isotropic GC systems; we use model GT for NGC 1407 
(see §4.1.4|) rather than the consensus model C. In the 
case of M87, the RKOl model has been used rather than 
the mass model of C-l-01 — which attempts to fit the over- 
all Virgo cluster dynamics but does not match well the 
steep profile increase in the region constrained by the X- 
rays and GCs. We do not include NGC 3379 and M60 
since clear results are not yet available based on the GCs 
alone^° . 

The summary of these results is shown in Fig. 1261 
In most of these systems, an upturn in Vc(r) is visi- 
ble: this is an intrinsic feature of ACDM halos that is 
masked in the central parts by the mass contribution of 
the galaxy, where it co nspires to produce a roughly con- 
stant Vcjr) profile (q.v. iGerhard et al.ll200]l iRusin et al.l 



[2003t iGavazzi et al.l 120071 ^. In most cases, the upturn 
is mild at the radii probed by the GC data, but in the 

20 For NGC 3379, the GCs are compatible with the PN- 
based result that the vdr) profile declines rapidly with radius 
(|BerKond et al. 2006; Douglas ct al. 2007), while for M60, the X- 
ray results imply a roughly constant Vc{r) profile to ~ 30 kpc (see 
H-l-08). 
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Fig. 26. — Mass modeling results for group/cluster-central galax- 
ies, based on GC dynamics. For clarity, no unc erta inties are shown. 
The linestyles correspond to galaxies as in Fig. 1231 The thicker seg- 
ments of the curves with terminating points show the approximate 
radial range constrained by the data. 

cases of NGC 1407 and M87, a strong upturn is seen. 
Note that the vb,ms{Rp) profiles (cf. Fig. [23]) are not 
straightforward indicators of the Vc(r) profiles, which can 
be inferred only with a dynamical modeling process that 
includes the slopes of the tracer density profiles and (ide- 
ally) the tracer orbit anisotropies. 

Before comparing NGC 1407 and M87 further, we note 
that their best-fit GC-based mass profiles are not unique, 
and in fact the two DM halo profiles are nearly consis- 
tent within the uncertainties. However, the independent 
constraints from X-rays (for M87) and group galaxies 
(for NGC 1407) do support these best-fit models. The 
difference in Vc amplitudes between the two systems im- 
plies a factor of 2 difference in DM content within ~ 30 
kpc (consistent with the observation that the X-ray tem- 
peratures differ by a factor of 2). The M87 sub-cluster 
appears to have a much higher halo concentration than 
the NGC 1407 group, which suggests that the NGC 1407 
group collapsed at a much lower redshift than than core 
of the Virgo cluster did. 

To first order, a CDM halo's initial collapse redshift Zc 
is related to its characteristic density Ps (ta bulated for 
NGC 1407 in Table [1 by oc (l-f z^)^ (Na varro et all 
119971: iBullock e t aT 2001). In more detail, the ACDM 
toy model of [Bullock ct al. (2001) implies Zc ~ 1.7 for 
NGC 1407 (age ~ 10 Gyr) and Zc ~ 4 for M87, assuming 
no alteration of the central DM density through baryonic 
effects^i. The halo mass profile of NGC 1399 is still fairly 
uncertain (see R-l-08), but the best-guess central mass 
profile implies Zc ~ 4. These numbers also reflect our re- 
sult of a "normal" DM halo concentration for NGC 1407, 
while the other systems have fairly high concentrations, 



2^ ITrentham et al] II2006I ) used different criteria to suggest 
Zc ~ 0.8 for the NGC 140 7 grou p, while the simulations 
of Ivon Benda-Beckmann et al.l 1)20081 ) predict Zc ^ —0.1 for 
NGC 1407, and z ~ 0.05 for M87 and NGC 1399. 
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since at fixed mass, higher-concentration halos will col- 
lapse first. (Note that if the X-ray result for NGC 1407 
is correct, then it has Zc ^ 5.) Optically, the Fornax 
core seems somewhat more relaxed than Virgo, with a 
higher galaxy density (also so mewhat higher than that 
of the NGC 1407 group; B-l-06b: , Jordan et al.ll2007t l. sug- 
gesting an earlier collapse — but the DM density should 
provide the more accurate picture. Thus the NGC 1407 
group appears to be at an earlier evolutionary stage than 
both Virgo and Fornax (i.e., starting its collapse later). 

5.3. Baryon content 

We next consider some global baryonic properties of 
the NGC 1407 group, comparing them to empirical and 
theoretical values for large samples of groups in the uni- 
verse. We start with the average overall mass-to-light ra- 
tio T for galaxy groups. There are many possible sources 
for these values, most of th em wit h similar findings, and 
we select the results of Eke et all ([2006), who used two 
different techniques to derive group and cluster T values 
based on - 100,000 2dFGRS galaxies at z - 0.1 (their 
halo mass definitions are within a few percent of ours 
for r^Af = 0.3). The first is a direct empirical technique 
of using the measured redshifts to calculate the masses 
of the individual galaxy groups with a virial analysis. 
The second is a hybrid empirical-theoretical matching of 
the observed luminosity function of groups to the mass 
function expected from ACDM (we are not aware of any 
purely theoretical prediction, owing to the highly uncer- 
tain effects of baryonic physics in structure formation). 
As displayed in Fig. 1271 both approaches yield very simi- 
lar results, with T increasing from a value of ~ 60 ^Q.bj 
for group luminosities of L*, to a maximum value of 
~ 300 T0^bj for clusters^^, which is similar to the av- 
erage T for the whole universe. 

We plot the results for the NGC 1407 group, increas- 
ing the B-band lumino sity by 10% to m ap to the pho- 
tographic 6j-band (see lEke et al.|[2004b[ ). Here we see 
that at the optical luminosity of the NGC 1407 group, a 
typical T ought to be ^ 130-190 Tq^;,^. Our consensus 
solution is much higher than this, with T ~ 800 Tq^^^, 
which is extreme even for a galaxy cluster. However, the 
range of uncertainty for NGC 1407 is still large, and we 
note that the shaded prediction region reflects the un- 
certainty in the mean value, not the cosmological scatter 
which may be larger. The NGC 1407 T values could be 
brought down if more luminosity were included, e.g., in 
additional group member galaxies, extended stellar ha- 
los, or intra-group light. However, any additional mem- 
bers will be dwarfs and not contribute much to the total 
luminosity, while the latter two issues should affect the 
2dFGRS survey results similarly. Note that in the B- 



and X-bands, the NGC 1407 group T is 850 
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and 1901|q° '^q,k^ respectively (for solution GG, where 
no X-ray constraints or concentration priors are applied, 
T is ~ 60% higher). 

The high T of the NGC 1407 group goes along with 
other unusual properties. The low luminosity for its halo 
mass means that it is also deficient in bright galaxies, at 
a level ~ 2-3 times lower than found observ ationally for 
typical low-redshift groups of similar mass (|Yang et al.l 



IMountrichas &: Shanks! II2007I ) used a difTerent analysis on the 
same data and found T higher by a factor of ~5. 
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Fig. 27. — Virial mass-to-light ratio, as a function of group 
luminosity. The dashed line with shaded region shows the empir- 
ical results (with uncertainties) from Eke et al. (2006) based on 
group dynamics. The solid line shows the semi-empirical result 
from the same authors, based on luminosity-mass function match- 
ing. The horizontal dotted line shows the mean value for the uni- 
verse, assuming a matter density of Cm = 0.3; a Hubble constant of 
Hq=70 km Mpc~^ is adopted. The solid points with error bars 
show results for the NGC 140 7 group (model C fr om this paper) 
and for the Fornax cluster tPrinkwater et 311120011 : R-l-08), based 
on models of their internal dynamics. The open points show more 
uniform results for the same systems and for the Virgo cluster, i.e., 
using the same s urvey data and mass analyses (from the projected 
mass-estimator; ITullvl [2005i). The star symbol shows a rough es- 
timate for the eventual Eridanus supergroup/cluster, extrapolated 
from B-|-06b. 



[2OO5I: [Hansen et all [20091 ) . although NGC 1407 appears 
to be an outlier at only the ^ 1-2 a level. The luminos- 
ity of the central galaxy NGC 1407 i tself is a lso low at 
the ^ 2 a le vel for groups of its mass (|Lin fc M ohr 2004 



lYang et al ' 2005: Brough et all 120081: 
but see IConrov fc Wechsleij |2009[ ) , 
type galaxy fraction of 



More et all 12009 



The group's early- 
80%-90% (B-h06a; B-K06b) is 
higher than typical for a group of its mass or velocity dis- 
persion (Zabludoff & Mulchaey 1998; Y ang et all l2005t 
B+06a). The only thing "normal" about the NGC 1407 
group is the luminosity of th e central galaxy re lative to 
the total luminosit y (OP04: [LhTfc Mohi|[200l B-h06a; 
[Hansen et alll2009[ ). 

For comparison, we also plot in Fig. [53 some literature 
results for the two nearest galaxy clusters, Fornax and 
Virgo. These suggest that Fornax is more than twice 
as luminous as the NGC 1407 group, yet ~ 30% less 
massive, resulting in a factor of 3 difference in T. Virgo 
on the other hand may have an T consistent with that of 
the NGC 1407 group, suggesting that these two systems 
are physically similar despite their very different galaxy 
content. 

One key uncertainty in these conclusions is the dis- 
tance: although this is ideally the era of distance mea- 
surements to 10% accuracy, NGC 1407 is a somewhat 
problematic case. As summarized by S-l-08, distance es- 
timates to this galaxy using various well-calibrated tech- 
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niques range between 19 and 27 Mpc. We have adopted 
21 Mpc since larger distances would seemin gly introduce 
proble ms for the GC luminosities and sizes (| Forbes et al.l 
l2006bf) . and incidentally would make it even harder to 
fit the X-ray data unless an implausibly low stellar T is 
adopted (see §4.3.ip . Most of our results in this paper 
have been independent of distance, e.g., the shape of the 
Vc{r) profile, but changing from our default distance of 
21 Mpc to 27 Mpc would increase the halo mass and 
shift T downward (though still leaving NGC 1407 as an 
outlier in Fig. [27)1 . 

An alternative way to consider the NGC 1407 group 
is via its baryon fraction /b = Mb/Mvir, where Mb 
is the total baryon mass including stars (in and be- 
tween galaxies), gas, and any other undefined bary- 
onic form. Assigning a conservative uniform value of 
T*,s — 4.45T0_B to all of the group galaxies, the to- 
tal stellar mass is ^ 4 xlO^^M©; the hot gas mass is 
~ 10-^°Mq (Z-I-07). Our consensus mass result then 
implies that .f t, ^ 0.002-0.007. G iven the cosmologi- 
cal /b = 0.17 (|Hinshaw et al.ll2009[) . this means '^96%- 
99% of the baryons in the system have been lost or 
else are invisible. Again, incorporating higher val- 
ues and "missed" galaxy halo light could boost /b, but 
only by a factor of 2. Such a result would contra- 
dict theoretical expectations for cl uster-mass halos to re- 
main essenti ally c l osed-boxes (e.g..lpe Lucia et al.]|2004l : 
lEttori et all l2006t ICrain et all 12007D, and observational 
evidence for bar yon "disappearance" t o be relatively con- 
stant at ~ 25% (I Gonzalez et al.ll2007[) . 

The under-represent ation of baryons may occur in 
many stellar systems (jMcGaughl |2008() . which might 
be explained by expulsion from feedback mechanisms — 
even in systems as massive as NGC 1407 (e.g., 
iFinoguenov et ai]|2001l ). Intra-group stellar light could 
also be considered a candidate for the missing baryons, 
except that this would comprise ~ 40 times as much stel- 
lar light as the rest of the group galaxies combined. A 
third alternative is that most of the missing baryons re- 
side i n a relati vely cool and/or diffuse gaseous medium 
(e.g. JCen fcO strikcr 2006; Nicastro ct al. 200_8.), which 
upon collapse of the supergroup into a cluster might heat 
up to emission levels detectable by X-ray telescopes. 

5.4. Implications 

Bringing together all of the strands from the preced- 
ing subsections, we consider their possible implications 
for the evolutionary history of NGC 1407 and its group, 
and where this system fits in a cosmological context. Re- 
visiting first the membership status of NGC 1400 with 
its large peculiar velocity, we note that there is no evi- 
dence from photometry and kinematics that it is associ- 
ated with a second group of galaxies superimposed along 
the line-of-sight. Also, the high inferred group mass does 
not hinge upon this single velocity, and exceeds the mini- 
mum value ( ^ 3 X lO^^Af©) necessary to keep NGC 1400 
bound. Interestingly, if NGC 1400 is now at pericenter 
on a near-radial orbit, then our best-fit solution "GT" 
implies that its apocenter is at rvir, and thus that this 
galaxy has only recently joined the group (as also sug- 
gested by S-l-08). 

Without this new arrival, NGC 1407 would techni- 
cally meet the definition for a fossil group, i.e. a sys- 
tem that consists optically of a single bright galaxy 



but has a group-mass halo (e.g., iJones et all l2003f) . 
The implication is that the GCG swallowed up any 
other bright primordial members at early times, re- 
sulting in a transient fossil phase as experienced by 
many grou ps (D'Onghia et al. 20051: iDariush et al.ll2007l : 
Ivon Benda-Be ckm ann et_ _al. .200^ . 

The GC kinematics in NGC 1407 and other GCGs 
(i.e., the moderate amount of rotation and flat veloc- 
ity dispersion profiles) appear broadly consistent with 
the formation of these galaxies in multiple major merg- 
ers of gas- free galaxies with pre-existing GCSs (SJSH]). 
However, the involvement of the group environment is 
unclear, and more observational and theoretical work 
is needed to determine, e.g., when the GCGs formed 
relative to their surrounding groups, and whether the 
GCSs are more closely associated w i th the galaxies or 
the groups ( d. iBlakeslee et al.l[l997L lMcLaughlii3[l99l 
iSpitler et al.ll2008bH . 

The detailed mass profiles derived from the GC kine- 
matics generally show an onset of the group- or cluster- 
scale halo via an increasing Vc{r). A comparison of the 
DM densities around NGC 1407, M87, and NGC 1399 
implies that the cores of the Virgo and Fornax clusters 
collapsed at earlier times than the NGC 1407 group did. 
Some commonalities in the GC kinematics and mass pro- 
files of NGC 1407 and M87 suggest that the formational 
pathways of these two systems were similar, albeit at 
different redshifts. The properties of the DM itself could 
in principle be further probed via the central slopes of 
the mass profiles and from deriving a mass-concentration 
relation — but any such conclusions will require greater 
care in dealing with selection effects, and with system- 
atic issues such as orbital anisotropy. 

Comparing the NGC 1407 group's global properties to 
empirical results for large samples of groups, we find that 
its high T and correspondingly low star formation make 
it a "dark cluster" as proposed by iGouldl (|1993( ) . The 
overall impression is of a poor galaxy group trapped in a 
high-mass halo that was somehow unable to form more 
galaxies, and which evolved its galaxies quickly to the 
red sequence owing to interactions or to the lack of fresh 
gas for star formation. The system could represent the 
tail end of galaxy cluster formation in the universe (with 
its peak at z ^ 0.7 for mas ses of - lO^^M^: [lacker et all 
120061: IConrov fc Wechsle"^l2009[ ). The future collapse of 
the Eridanus "supergroup" into a cluster should bring the 
system's T closer to typical values (see Fig. [27| . while the 
missing baryons might emerge as the hypothetical unseen 
gas is heated into X-ray bright gas. The only incongruity 
would be the unusually low luminosity of NGC 1407 as a 
cluster-central galaxy, unless it experiences another ma- 
jor merger. 

One could posit that high-T clusters such as Virgo 
were built up from progenitor groups like NGC 1407 with 
high T (but earlier collapse times), akin to the notion 
that m any brightest cluster gala xies originated as fossil 
groups (jKhosroshahi et all 12006). This scenario has in- 
triguing relevance for the observational puzzle that the 
stellar mass fraction in low-redshift galaxy groups and 
clusters decreases with mass — contradicting the expecta- 
tion from CDM hierarchical merging that more massive 
systems a re formed from less massive systems with simi- 
lar (.Balogh et al., ,2008). The NGC 1407 group on the 
other hand has a very low ^ 0.004, which would go 
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in the direction needed by theory. An obvious interpre- 
tation is that this group is a reUc from an earlier epoch, 
when the progenitor groups of clust ers had systernaticall y 
lower However, as discussed bv lBalogh et al.l ()2008D . 
the other present-day groups should also have low /, un- 
less they had problematically recent bouts of star forma- 
tion. One way out is to suppose that NGC 1407 rep- 
resents the tip of an iceberg, of groups that were either 
missed observationally, or had their masses systemati- 
cally underestimated. A less mundane alternative is that 
structure formation on these scales is anfi-hierarchical 
(e.g., because of warm DM). 

The notion that clusters are built up from a biased 
subset of groups would also point to a pre-processing 
scenario for cluster galaxy evolution. However, the trend 
of recent findings is that galaxy clusters have fairly uni- 
versal properties related to their total mass rather than 
to their merger histories, i.e. the galaxy content is some- 
how systematically altered as the clusters are assembled, 
rather than the clusters be ing as sembled from sys temat- 
ically distinct groups (e.g ., lLin et al. 2006; M uzzin et al.l 
I2007t [Hansen et all [20091) . Add itional nuances are be- 
ing added to this pic ture (e.g., IWeinmann et al.l 120061 : 
iPoggianti et al.l [2006f ). and it is possible that a mass- 
correlated parameter such as density could be the main 
driving factor. In this context, the NGC 1407 group 
could be a telling system, since it is an outlier in the usual 
optical correlations, and its mass is bel ow the proposed 
boun dary for galaxy quenching to occur ([Poggianti et all 
[2006[) . 

To better understand the implications of the NGC 1407 
group, its uniqueness must be clarified. There are obser- 
vational problem s with detecting such systems, as dis- 
cussed bv lGouldl ([l993h : the "finger of God" effect of a 
high-mass poor group tends to be lost in the crow d of 
neighboring redshifts (see also Cow an fc Ivezidl2008l ). It 
would thus not be surprising if many such systems have 
been missed in optical surveys. On the other hand, the 
group is easily detected in X-ray surveys, with a rough 
Tx-based guess for T that is only a factor of ~ 2 lower 
than our detailed estimate (see §4.2|) . Furthermore, it 
is an ordinary system when compared to X-ray corre- 
lations involving Lx, T x, Cp, pspec, Px, M^-„, Lb, etc. 
([O'Sullivan et al.ll200H OP04; B-h06a). Its only unusual 
aspect seems to be a Tx somewhat higher than typical 
for its Lx, although the Tx(r) profile in detail may not 
be unusual for a group asmussen fc Ponmanl[2?)07[ : cf. 
iDiehl &: Statleil[2008af ). This contrasts with its optical 
properties, almost all of which appear to be unusual (see 
§5.3p . There may be something of a contradiction be- 
tween optical- and X-ray-based mass functions which we 
cannot resolve in the present paper (cf. iMendel et al.l 
[20081) . 

There is anecdotal evidence for systems similar to the 
NGC 14 07 group. T he fairly small nearby group cat- 
alogs of iTullvl (|20050 and B+06a include a few more 
candidates for ultra-dark groups — suggesting that the T 
of NGC 1407 may be well within the cosmological scat- 
ter. Two low er-mass groups at inte rmediate redshift in 
the survey of I Poggianti et ahl (|2006( ) were found to have 
unusually high early-type galaxy fractions, prompting 
the authors to dub them "bare" massive-cl uster cores. 
The "dark core" in the cluster Abell 520 (Mahdavi et al.l 
|2007( ) could be a preserved descendant of a group like 



NGC 1407, although this core appears to be missing its 
GCG. High-T systems could also be detected as "dark 
clumps" in weak lensing su rveys (cf. lErben et al.ll2000l : 
Ivon der Linden et al.ll2006f) . 

Another fairly nearby group associated with NG C 1600 
may a lso have Tb ^ 1000 T^,©, as found bv Smith et al.l 
(|2008D . who suggest that this system represents a new 
class of fossil group with relatively low Lx- Alterna- 
tively, NGC 1407 and NGC 1600 might be heralds of a 
lurking population of massive ultra-dark groups, which 
were recognized in these cases only because of their bright 
GCGs and X-ray emission. A careful mining of modern 
observational surveys, and of cosmological simulations, 
may reveal the frequency of such systems, and their place 
within the wider context of structure formation. 

6. CONCLUSIONS 

We have used Keck/DEIMOS to study the kinematics 
of GCs around NGC 1407, arriving at a final sample 
(including literature data) of 172 GC velocities out to 
galactocentric radii of ~ 60 kpc. Outside the central 
regions, the GC system has weak rotation (urot/o'p ^ 0.2) 
around the major axis, with the metal-poor and metal- 
rich GC subsystems showing rotational misalignments. 
We also detect a moving group of GCs at large radii. 

The GC velocity dispersion profile declines rapidly to 
~ 2', whereafter it stays constant or rises gradually. The 
rapid decline may be due to a contribution of centrally- 
located unresolved dwarf- globular transition objects on 
high-eccentricity orbits; we find two clear photometric 
examples of such objects. We also identify one probable 
intra-group GC, and find that the "blue tilt" (a color- 
magnitude relation for metal-poor GCs) discovered with 
HST imaging persists with ground-based photometry of 
a spectroscopically-confirmed GC sample. 

Dynamical modeling of the GCS indicates a circular 
velocity profile that continues to rise out to at least 
60 kpc, consistent with a massive ACDM halo with 
Afvir ^ 7 X IQ^^Mq and a group mass-to-light ratio of 
Tb ^ 800 T0^B. This agrees with results from group 
galaxy kinematics, but our preliminary X-ray based anal- 
ysis of the mass profile shows some discrepancies. In 
order to reach agreement, the GCs would have to fol- 
low an implausible orbital distribution, with highly ra- 
dial anisotropy in the center rapidly transitioning to tan- 
gential anisotropy in the outer parts. We survey X- 
ray results in other systems, finding some suggestions 
of discrepancies with kinematics-based results. Also, 
the common phenomenon of X-ray "kinks" implies prob- 
lems either with the assumption of hydrostatic equilib- 
rium, with standard stellar T values, or with the CDM 
paradigm. Resolving this issue is of obvious importance 
to cosmological inferences from X-ray-based mass func- 
tions. 

We compare NGC 1407 and other nearby group- and 
cluster-central galaxies, finding similarities in their GC 
kinematics and mass profiles. We suggest that the 
NGC 1407 group is similar to, but collapsed later than, 
the "seed" progenitor group that became the core of the 
Virgo cluster. Comparing its properties to large surveys 
of galaxy groups and clusters, this group appears to be a 
strong outlier. It is one of the "darkest" systems in the 
universe, whose T is simi l ar to those of the faintest dwarf 
galaxies ([Simon fc Gehal[2007[) . and we are left with a 
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puzzle as to why star formation was so attenuated in 
this system, and where the initial cosmological baryon 
complement is hiding. Finding more such high-T groups 
is part of the difficult but crucial inventory of mass and 
light in the universe, and could help disentangle the roles 
of halo density and mass in driving galaxy evolution. 

Our new observations compound the existing evidence 
that GC systems have isotropic or tangentially-biased 
velocity distributions, which seems to dramatically con- 
tradict most theoretical expectations for halo particles 
to reside on radially-based orbits. However, there are 
open questions about the orbital dynamics of GC acqui- 
sition and destruction in a realistic theoretical context, 
and about possible observational links between GC lu- 
minosities, sizes, and orbits. 

NGC 1407 and its surrounding galaxy group are an 
intriguing system that warrants further investigation. 
More complete kinematics observations and more de- 
tailed dynamical models are necessary to unravel its or- 
bital and mass-to-light ratio mysteries. Our efforts con- 
tinue along these lines. 
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